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Abstract 

A  survey  Is  made  of  various  singular  and  nonuniform 
limits  of  solutions  of  the  Boltzmann  equation.   Included 
are  the  limits  of  large  and  small  mean  free  path,  fast 
and  slow  molecules.  Initial,  boundary,  and  shock  layers, 
and  singular  features  of  the  collision  Integral.   Several 
new  results  are  given  Including  an  analysis  of  the 
Hllbert  expansion  In  steady  flow,  of  the  kinetic  boundary 
layer  Including  a  complete  formal  theory  of  slip  coeffi- 
cients, and  of  the  profile  of  a  very  strong  shock. 
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1 .   Introduction 

The  study  of  singularities  and  nonuniform  limits 
appears  to  be  a  major  preoccupation  of  applied  mathemati- 
cians.  Aside  from  the  intrinsic  Interest  of  such  problems, 
there  is  good  reason  for  this  concentration.   A  good 
existence  theorem  will  include  as  a  corollary  continuous 
or  smooth  dependence  of  the  solution  on  parameters.   If  a 
parameter  which  occurs  in  the  problem  has  a  natural  range 
of  variation,  then  one  can  expect  the  behavior  at  the  end- 
points  to  be  singular;  otherwise  the  solution  could  be 
continued.   Thus  a  detailed  description  of  the  singular 
endpoints, together  with  continuous  dependence  in  between, 
should  give  a  uniformly  valid  qualitative  picture  of  the 
entire  problem. 

There  is  also  a  quantitative  profit  that  can  be  gained 
from  knowledge  of  the  singular  endpoints.   It  is  frequently 
possible  to  "subtract  off"  the  singular  part  and  then 
obtain  quite  accurate  answers  for  the  remainder  using 
relatively  unsophisticated  or  ad  hoc  procedures.   In 
kinetic  theory,  this  combination  of  relatively  recent 
precise  investigations  of  many  special  limiting  cases  with 
the  older  moment  and  polynomial  approximations  offers  the 
likllhood  of  quite  accurate  practical  calculations  even  if 
convergence  is  not  yet  established. 
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In  a  previous  talk  at  this  Symposium  ,  existence 
theorems  covering  a  certain  area  of  kinetic  theory  were 
given  together  with  a  rigorous  analysis  of  the  behavior 
at  one  singular  limit  (small  mean  free  path).   We  now 
survey  all  (or  at  least  many)  of  the  singular  limits  to 
which  solutions  of  the  Boltzmann  equation  fall  prey. 
Since  the  available  existence  theory  does  not  cover  this 
extended  territory,  the  present  treatment  will  be  formal 
rather  than  rigorous. 

For  coherence  in  the  oral  presentation,  the  emphasis 
was  placed  on  a  broad  survey  rather  than  on  the  new  material 
We  reverse  this  emphasis  in  the  written  text  since  there 
are  other  sources  available  for  the  details  of  many  of  the 
nonuniform  limits  which  are  here  coordinated.    The  new 
material  consists  of  a  unified  formal  treatment  of  the 
small  mean  free  path  limit  in  steady  boundary  value  problems 
(Sees.  7-10)  and  an  analysis  of  the  profile  of  a  very  strong 
shock  (Sec.  11). 

The  boundary  value  problem  separates  into  an  analysis 
of  the  interior  and  an  evaluation  of  asymptotic  boundary 
conditions.   For  the  interior  analysis,  an  asymptotic 


"Asymptotic  Eauivalence  of  the  Navier-Stokes  and  Non- 
linear Boltzmann  Equations,"  New  York,  April  22,  1964;  pub- 
lished iin  Proceedings  of  Symposia  on  Applied  Mathematics 
(American  Mathematical  Society,  I965),  Vol.  XVII. 

References  will  be  given  as  appropriate. 
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procedure  Is  given  which  is  superior  to  the  direct  Hllbert 
expansion  or  the  conventional  Chapman -Enskog  expansion; 
(but  In  special  cases  this  reduces  to  a  degenerate  form 
of  the  Hllbert  expansion).   This  analysis  of  the  Interior 
of  the  domain  Is  supplemented  by  a  formally  complete  des- 
cription of  the  slip  flow  boundary  conditions  In  a  general 
domain.   Numerical  evaluation  of  the  slip  coefficients 
requires  solution  of  certain  special  one -dimensional  prob- 
lems generalizing  the  Kramers  problem. 

The  strong  (strictly  speaking.  Infinite  strength)  shock 
profile  Is  reduced  to  a  relatively  nonslngular  problem  by 
Introducing  the  correct  scaling  In  physical  space  and  In 
velocity  space  after  subtracting  off  a  suitable  singularity. 
This  treatment  can  be  qualitatively  Interpreted  as  a  more 
rational  refinement  of  the  ad  hoc  Mott -Smith  blmodal 
approximation.   In  the  correct  scaling.  It  appears  that  a 
very  simple  Interpolation  can  be  made  between  the  Infinitely 
strong  and  weak  shock  profiles. 

A  survey  of  various  singular  limiting  cases  Is  pre- 
sented In  Sees.  2-6,  followed  by  more  detailed  consideration 
of  steady  slip  flow  In  Sees.  7-10  and  of  the  shock  layer  In 
Sec.  11.   We  conclude  In  Sec.  12  with  some  opinions  on 
future  trends  In  kinetic  theory. 

We  should  like  to  acknowledge  many  useful  discussions 
with  Dr.  Yoshlo  Sone  In  connection  with  the  material  in 
Sec.  10. 
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2.   Nonuniform  Limits 

The  argument  of  the  Boltzmann  equation  Is  a  distribu- 
tion function,  f{i,x,z).      This  represents  the  density  at  a 
given  Instant  In  the  phase  space  of  a  molecule;  x  =  position 
and  ^  =  velocity.   We  find  nonuniform  limits  at  each  extreme 
of  every  argument  of  f .   There  are  nonuniform  limits  as 
X  — >  0  (boundary  layer)  and  as  x  — >  oo;  also  as  t  -^   0 
(Initial  layer)  and  as  t  — >  oo  (approach  to  equilibrium); 
also  for  slow  molecules,  i   — >  0,  and  fast  molecules,  ^  — >  «. 
If  dlmenslonless  variables  are  Introduced  for  i,   x,  and  t 
a   relative  to  a  representative  mean  thermal  speed,  x  and  t 
relative  to  some  given  scale)  ,  then  a  single  dlmenslonless 
parameter,  e,  appears  In  the  Boltzmann  equation.   This  Is 
the  ratio  of  mean  free  path  to  the  given  scale  length  or, 
equlvalently,  of  the  mean  collision  time  to  the  scale  time. 
Both  limits,  e  — >  0  (continuum)  and  e  — >  oo  (free  flow)  are 
singular. 

Since  macroscopic  fluid  dynamics  occurs  as  a  limiting 
case,  e  — >  0,  In  the  theory  of  the  Boltzmann  equation,  all 
singular  and  nonuniform  limits  of  classical  fluid  dynamics 
are  automatically  Included.  These  encompass  an  assortment 
of  boundary  layers,  wakes,  shock  layers,  various  lineariza- 
tions, behavior  at  Infinity,  small  and  large  Mach  number 


We  prevent  the  appearance  of  a  Mach  number  In  the 
equation  by  taking  the  i   scale  equal  to  the  ratio  of  the 
X  and  t  scales. 
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and  Reynolds  number,  and  many  more.   The  singular  approach 
to  fluid  dynamics  as  e  — >  0  is,  of  course,  supplemental 
to  the  strictly  fluid  nonunlf ormltles . 

An  Interesting  point  Is  the  fact  that  the  parameter  e 
which  gives  rise  to  fluid  dynamics  In  the  limit  e  — >  0 
occurs  within  fluid  dynamics  proper  as  a  finite  parameter, 
e  =  M/Re  (Knudsen  number  equals  Mach  over  Reynolds  number). 
Prom  fluid  dynamics  alone  we  could  not  suspect  that  finite 
(i.e.,  not  small)  values  of  I^Re  are  Intrinsically 
unreliable.   The  limit  e  — >  0  Is  singular  In  both  fluid 
dynamics  and  kinetic  theory,  with  the  two  sets  of  compli- 
cations superposed.   For  example,  a  fluid  dynamic  boundary 
layer  expansion  In  e  (Re  — >  oo  at  finite  M)  requires 
supplementation  by  kinetic  effects,  usually  at  second 
order  In  e . 

The  fact  that  fluid  dynamics  emerges  from  the  Boltzmann 
equation  In  a  certain  limit  Is  Intimately  tied  to  the  laws 
of  conservation  of  mass,  momentum,  and  energy.   The  con- 
servation laws  result  from  certain  elementary  properties 
of  the  collision  Integral,  viz.  that  the  nullspace  is 
nontrivlal  (and  finite  dimensional).   This  situation  can 
be  compared  to  the  case  of  neutron  transport  where  there 

■X- 

is  a  similar  equation,  also  called  the  Boltzmann  equation. 


*  „ 

Cf.  N.  Corngold,   Some  Recent  Results  in  the  Theory 

of  the  Transport  of  Thermal  Neutrons,"  to  appear  in 

Symposia  in  Applied  Mathematics  (American  Mathematical 

Society),  Vol.  XX. 
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In  the  neutron  case,  conservation  does  not  hold.   As  a 
result   solutions  will.  In  general,  either  grow  exponen- 
tially or  decay  exponentially  on  the  time  scale  of  the 
mean  collision  time.   To  identify  the  borderline  case 
between  growth  and  decay  (e.g.,  in  its  dependence  on 
geometry)  is  one  of  the  basic  problems  in  neutron  trans- 
port.  In  a  sense,  the  Boltzmann  equation  for  a  gas  (or 
Boltzmann' s  equation)  represents  a  set  of  measure  zero 
from  the  point  of  view  of  neutron  transport.   But  it  is 
exactly  this  borderline  quality  which  gives  a  gas  its 
unique  and  distinctive  features  (including  fluid  flow) 
and  makes  a  significant  portion  of  the  mathematical  analysis 
much  more  delicate. 

Our  primary  concern  is  with  nonunif ormities  which 
are  peculiar  to  kinetic  theory  (one  of  which  is  the  singu- 
lar connection  between  the  fluid  and  kinetic  descriptions, 
cf.  Sec.  5).   Many  nonunif ormities  arise  from  the  specific 
form  of  the  kernel  of  the  collision  term  of  the  Boltzmann 
equation  which  represents  the  nature  of  the  intermolecular 
force.   It  is  Important  to  realize  that  "the"  Boltzmann 
equation  is  not  a  specific  equation  but  a  generic  term  for 
a  class  of  equations  (similar  in  extent  to  the  term  integral 
equation) .   There  are  various  types  of  collision  kernels, 
some  more  singular  than  others.   In  particular,  one  feature 
that  can  be  extracted  from  the  collision  kernel  is  the 
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velocity-dependent  collision  frequency  v(^)  and  the 
related  velocity-dependent  mean  free  path  A(^)  =  ^/v(4). 
Much  of  the  singular  behavior  for  small  and  large  ^  can 
be  related  to  the  corresponding  limiting  behavior  of 
v(^)  and  A(^),  cf.  Sees.  3  and  ^. 

One  class  of  collision  kernels  about  which  almost 
nothing  is  known  (in  a  basic  sense  the  most  important 
class)  is  that  for  a  non-cutoff  Intermolecular  potential. 
In  this  case  no  simple  collision  frequency  v(^)  or  mean 
free  path  A(^)  can  be  defined,  and  the  nature  of  the  col- 
lision operator  (e.g.,  whether  it  is  an  Integral  operator, 
singular  integral  operator,  differential  operator,  etc.)  is 
not  known.    For  a  certain  type  of  cutoff  (in  angle  of 
deflection),  a  relatively  complete  mathematical  theory  has 
been  developed  [Grad,  1963b].   For  another  cutoff  (inter- 

jt  jt  jt 

molecular  distance)  a  few  results  have  been  obtained. 
One  would  like  to  know,  on  physical  grounds,  that  "many" 
properties  of  solutions  of  the  Boltzmann  equation  are 
insensitive  to  the  type  of  cutoff  and  are,  indeed,  similar 


For  a  result  in  a  special  case  see  [Smolderen,  I965; 
Willis,  1965]. 

An  approach  to  this  problem  is  described  in 
[Finkelstein,  I965]. 

*** 

C.  Cercignani,  private  communication. 
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to  those  for  non-cutoff  potentials.   The  one  case  In  which 
this  Is  known  to  be  true  Is  In  the  formula  for  the  calcu- 
lation of  transport  coefficients  [Burnett,  I956].   What 
is  probably  the  case  more  generally  Is  that  the  precise 
nature  of  certain  mathematical  singularities  Is  sensitively 
dependent  on  the  nature  of  the  cutoff  (or  on  the  exact 
decay  of  a  non-cutoff  potential),  even  though  some  gross 
properties  of  the  gas  flow  will  not  be  so  sensitive. 

We  have  mentioned  that  all  the  singular  limits  of 
classical  fluid  dynamics  are  contained  within  kinetic  theory. 
This  Is  not  quite  so;  some  of  the  bad  behavior  may  be  wiped 
out.   In  particular,  nonlinear  and  linearized  results  are 
related  quite  differently  In  the  fluid  and  kinetic  versions. 
The  streaming  term  of  the  Boltzmann  equation  Is  strictly 
linear;  nonllnearlty  enters  only  through  the  quadratic 
collision  term.   But  fluid  equations  are  obtained  from  the 
Boltzmann  equation  by  a  procedure  which  operates  only  on 
the  linear  streaming  terms,  explicitly  Ignoring  the  colli- 
sion term;  the  relevant  moments  of  the  collision  term 
vanish  because  of  conservation.   The  complete  set  of  fluid 
conservation  equations  Is  therefore  linear  (in  the  proper 
variables)  until  It  Is  "simplified"  by  use  of  the  strongly 
nonlinear  Maxwelllan  distribution  function.   The  nonllnearlty 
of  the  fluid  equations  (and  the  associated  equations  of 
state)  Is  Imposed  by  the  nonlinear  Boltzmann  collision 
term,  but  only  In  a  rather  devious  manner. 
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To  Illustrate  this  divergence  between  the  two  theories, 
consider  the  flow  at  finite  Mach  number  past  a  small  sphere 
or  cylinder  (small  compared  to  a  mean  free  path).   The 
fluid  solution  is  definitely  nonlinear  regardless  of  the 
size  of  the  body.   The  kinetic  solution  also  yields  a 
strong  disturbance  near  the  body.   But  this  strong  distur- 
bance is  governed  by  the  strictly  linear  streaming  terms 
of  the  Boltzmann  equation;  the  collision  term  is  negligible 
until  about  a  mean  free  path  away,  at  which  distance  the 
disturbance  has  become  small.   The  limiting  flow  is  given 
as  a  finite  amplitude  solution  of  a  linear  equation  (free 
flow).   The  leading  term  of  the  departure  from  free  flow 
is  nonlinear  in  three  dimensions  but  is  linear  in  two 
dimensions. 

A  similar  example  is  given  by  high  frequency  sound, 
e.g.,  as  generated  by  an  oscillating  wall.    The  conven- 
tional fluid  linearization  which  takes  derivatives  as  well 
as  fluid  perturbations  to  be  small  becomes  increasingly 
poor  as  the  frequency  and  wave  number  become  large.   But 
in  the  kinetic  analysis,  a  disturbance  with  frequency 
which  is  large  compared  to  the  collision  frequency  is 
dominated  by  the  streaming  terms,  and  the  nonlinear  collision 


*  r 

For  the  solution  with  a  relaxation  model  see  [Weltzner, 

1965]  and  [Mason,  I965];  with  the  Boltzmann  equation  at  high 
frequency  see  [Grad,  1966]. 
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term  offers  only  a  minor  correction  to  the  behavior  of 
even  a  finite  amplitude  wave.   For  a  fixed  small  ampli- 
tude of  the  disturbance,  linearization  of  the  Boltzmann 
equation  will  be  most  unreliable  at  some  intermediate 
frequency. 

A  similar  conclusion  is  reached  for  the  initial  value 
problem.   For  the  Boltzmann  equation,  the  natural  time 
scale  is  a  mean  collision  time.   One  would  expect  the 
discrepancy  between  a  linear  approximation  and  the  correct 
nonlinear  solution  to  grow  with  time  on  this  time  scale. 
But  if  the  disturbance  is  very  long  compared  to  the  mean 
free  path,  a  much  longer  (macroscopic)  time  scale  emerges. 
These  expectations  are  verified  by  quantitative  estimates 
in  linear  and  nonlinear  existence  theorems  [Grad,  1965]- 
Thus  the  worst  case  with  regard  to  linearization  is, 
again, the  intermediate  region  with  wave  length  comparable 
to  (or  probably  somewhat  larger  than)  the  mean  free  path. 

The  Navier-Stokes  level  of  the  Chapman -Enskog  expan- 
sion amounts  to  a  linearization  in  velocity  space  (small 
deviation  from  local  Maxwellian)  in  a  problem  which  may  be 
strongly  nonlinear  in  physical  space.   But  even  the 
linearization  in  velocity  space  is  not  uniformly  valid 
for  large  velocity  [Naraslmha,  I967]  (as  evidenced  by  the 
distribution  function  becoming  negative).   Even  In  a 
problem  in  which  the  fluid  disturbance  is  small,  i.e.. 
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linear  In  physical  space,  the  correct  treatment  of  the 
high  velocity  tall  Is  always  nonlinear;  the  tall  of  the 
distribution,  although  small,  Is  not  small  relative  to 
a  Maxwell Ian. 

Nonllnearlty  can  even  Introduce  a  dlsslpatlve  mechanism 
where  there  Is  none  linearly,  in  the  limit  of  small  mean 
free  path,  e,  one  aspect  of  dissipation  occurs  rapidly 
on  a  time  scale  e,  during  which  the  Boltzmann  H-f unction 
becomes  approximately  equal  to  the  thermodynamic  entropy  S; 
this  Is  then  followed  by  a  much  slower  viscous  dissipation 
on  a  time  scale  l/e-   The  linear  problem  Is  non-dlsslpatlve 
on  the  Intermediate  time  scale  0(1).   But  in  a  nonlinear 
problem,  shock  waves  provide  additional  dissipation  (which 
involves  H  and  S  differently)  on  the  time  scale  0(1 ). 

The  moment  and  polynomial  approximations  which  were 
originally  introduced  to  break  away  from  the  fluid  regime 
[Grad,  19^9]  are  essentially  ad  hoc.   They  usually  supply 
fairly  good  approximations  in  a  low  order  approximation, 
but  no  rapid  improvements  are  found  even  at  considerably 
extra  labor  [Pekeris,  et  al.,  I962].   But,  if  an  explicit 
singularity  or  nonuniformlty  which  prevents  rapid  conver- 
gence is  recognized,  the  polynomial  expansion  can  be 
altered  to  suit,  and,  although  still  ad  hoc ,  the  level 
of  accuracy  may  be  greatly  improved  even  at  low 
order.   The  most  elementary  example  is  the  inclusion  of 
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two-sldedness  In  the  distribution  by  using  half  range 
expansions  [Gross  and  Zlering,  1958;  Lees  1959]- 

Another  example  is  given  by  the  fundamental  solution 
in  steady  flow.   For  a  relaxation  model  of  the  Boltzmann 
equation  the  Fourier  transform  of  the  solution  is 
explicit  [Grad,  1959]-   For  the  actual  Boltzmann  equation 
we  are  not  so  fortunate.   However,  the  singularity  near 
the  source  i_s  explicit  to  a  certain  number  of  terms  in 
an  expansion.   Subtracting  this  off  leaves  a  residue  which 
is  bounded  and  can  be  expected  to  be  well -approximated  by 
a  simple  moment  expansion  (Euler,  Navier-Stokes ,  or 
Thirteen  Moment).   In  addition,  the  leading  terra  at 
infinity  is  governed  by  fluid  dynamics  and  is  therefore 
given  correctly  by  a  moment  approximation.   This  scheme 
has  been  carried  out  [Yu,  I967];  unfortunately  the  singu- 
larity at  the  origin,  although  "explicit",  involves  very 
complex  high-dimensional  integrals. 

A  more  successful  application  of  this  technique  is 
to  the  strong  shock.   We  recall  that  for  a  Navier-Stokes 
prof lle^  proper  scaling  (viz.  introducing  the  downstream 
mean  free  path  as  a  reference  length)  gives  a  shock  thick- 
ness which  does  not  vary  by  more  than  a  factor  two  or 
three  and  a  profile  which  is  not  very  far  from  universal 
as  the  strength  varies  from  zero  to  infinity  [Grad,  1952]. 
We  now  find  for  the  Boltzmann  equation,  after  subtracting 
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off  an  explicit  singularity  In  velocity  space,  that  there 
Is  a  definite  limiting  profile  in  ^  as  well  as  x  for  the 
distribution  function  at  infinite  strength.   This  limiting 
distribution  function  can  be  approximated  by  relatively 
crude  methods  (e.g.,  by  a  local  Maxwellian,  which  gives 
no  result  at  all  if  the  singularity  is  not  removed).   One 
conclusion  is  that  the  shock  profile  can  be  qualitatively 
approximated  over  the  entire  range  from  zero  to  infinite 
strength  by  simple  interpolation  (Sec.  11). 

This  very  brief  summary  of  singular  limits  is  con- 
tinued with  a  slightly  more  elaborate  discussion  of  the 
nonuniformlties  which  are  Introduced  by  the  collision 
kernel  and  at  the  large  mean  free  path  limit  in  the  next 
two  sections. 
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3-   The  Spectrum- -Fast  and  Slow  Molecules 

With  the  exception  of  the  finite  nullspace  of  the  col- 
lision operator  (which  expresses  conservation  and  leads  to 
fluid  dynamics),  probably  the  single  most  Important  quali- 
tative feature  of  the  linear  Boltzmann  equation  is  the 
almost  universal  occurrence  of  a  continuous  spectrum. 
There  are  two  sources  of  a  continuum,  the  streaming  opera- 
tor, ^'V^  and  the  collision  operator  L. 

The  streaming  operator  alone,  as  In 

has  the  entire  Imaginary  axis  as  a  continuum.   This  leads 
to  the  decay  In  time  of  an  arbitrary  Initial  function 
f(^,x,0)  [Grad,  I96I]  and  the  decay  In  space  of  the  signal 
induced  by  an  oscillating  wall  In  a  colllslonless  gas 
[Kahn  and  Mlntzer,  I965].   It  Is  also  the  dominant  effect 
at  high  frequency  even  In  the  presence  of  collisions. 
The  collision  operator  by  Itself, 

11+  L(f)  =  0  (5.2) 

has  a  continuum  covering  part  of  the  real  axis  for  most 
molecular  force  laws  which  have  been  Investigated  up  to 
the  present  time  [Grad,  1963b]. 

Since  the  operator  L  Is  symmetric  and  ^-V  Is  anti- 
symmetric, and  both  operators  are  unbounded,  the  spectrum 
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of  the  sum,  ^ -V  +  L,  Is  not  evident.   There  are  a  large 
number  of  Interesting  perturbation  problems  with  a  few 
proved  results  and  many  plausible  conjectures.   It  Is 
relatively  easy  to  make  reliable  plausible  estimates  of 
the  continuous  spectrum.   This  yields  valuable  and  detailed 
qualitative  Information  about  solutions  of  the  Boltzmann 
equation.   It  Is  also  the  most  Important  single  criterion 
for  assessing  what  phenomena  are  likely  to  be  misrepre- 
sented by  a  simplified  relaxation  model  of  the  Boltzmann 
equation.   In  particular,  many  singular  features  are 
Immediately  visible  from  a  study  of  the  continuous  spectrum. 

For  the  initial  value  problem  (3-2),  the  continuum 
Is  the  set  of  values  taken  by  the  collision  frequency  v(^). 
For  a  class  of  potentials  called  hard  we  have  v(0)  =  1 
Increasing  to  v(oo)  =  oo;  for  a  soft  potential,  v(o)  =  1 
decreases  to  v(oo)  =  0.   In  a  gas  of  hard  molecules  the 
high  velocity  tail  of  an  initial  distribution  f(^,0)  will 
Instantly  equilibrate,  this  effect  gradually  penetrating 
towards  the  center.   On  the  other  hand,  a  soft  gas  will 
reach  equilibrium  in  the  center  of  the  distribution  in  a 
time  0(1),  and  this  will  spread  more  and  more  slowly  to 
the  tail. 

A  more  interesting  illustration  is  given  by  a  boundary 
value  problem. 
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^1 1% -*- ^(^)  =  ^  (5-5) 


The  continuous  spectrum  (plausible  conjecture)  is  In  this 
case  determined  as  the  set  of  values  taken  by  the  recipro- 
cal of  the  free  path,  l/X^ii)  =  v  {i)/i-,-      In  most  cases  the 
continuum  covers  the  entire  real  axis.   Because  A,  =  0  for 
4-,  =  0,  there  Is  nonuniform  convergence  In  velocity  space 
for  the  distribution  function  near  a  wall.   The  difference 
between  f(4,x)  and  the  emitted  wall  distribution  f  {i ,0) 
(i     >  0)  does  not  approach  zero  In  a  velocity  boundary 
layer  of  thickness  I-,  =  0(x).   As  a  consequence,  certain 
moments  (e.g.,  temperature,  tangential  velocity)  have 

■X- 

Infinite  slopes  at  the  wall. 

At  the  opposite  extreme  we  have  A  — >  oo  for  i   — >  oo 
(with  the  single  exception  of  hard  spheres).   In  this  case, 
fast  particles  will  exert  an  Influence  at  long  distances. 
It  has  been  found,  for  example,  that  the  weak  shock  pro- 
file has  a  nonexponentlal  decay  at  Infinity,  apparently 
contradicting  the  Navler-Stokes  result  [Lyubarskl,  I961]. 
The  reason  Is  the  large  penetration  of  fast  particles 
past  the  bulk  of  the  shock  transition.   There  Is  an  Inver- 
sion of  limits  In  that  the  Navler-Stokes  profile  Is 


For  a  relaxation  model  see  [Sone,  196^a,b];  for  the 
Boltzmann  equation  see  [Grad,  I966]. 
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correct  for  a  relatively  greater  portion  of  the  profile 
as  the  shock  strength  decreases.   The  effect  Is  small 
because  there  are  exponentially  few  fast  particles. 
Whether  this  effect  is  important  in  a  strong  shock  has 
not  yet  been  determined. 

In  the  presence  of  a  wall  oscillating  at  a  fixed 
frequency, 

Icof  +  ^1  |i-+  L(f)  =  0  (3.^) 

the  spectrum  [i.e.,  of  Sf/^x-,  =  ^(f)]  moves  into  the 
complex  plane.   Many  qualitative  features  of  the  solutions 
can  be  predicted  just  from  a  study  of  the  continuous  spec- 
trum, e.g.,  the  existence  for  large  co  of  sub-boundary 
layers  in  the  approach  of  f  to  the  wall  distribution  f  , 
and  the  disappearance  of  a  sound  wave  (i.e.,  of  all  dis- 
crete normal  modes)  for  sufficiently  high  oo. 

The  situation  is  qualitatively  similar  (e.g.,  with 
regard  to  the  complex  continuous  spectrum  swallowing  up 
the  normal  modes  at  small  wavelength)  for  an  initial  value 
problem  at  fixed  wavelength, 

11+  le^k^f  +  L(f)  =  0  (3-5) 

A  more  detailed  description  of  these  spectra  is  dis- 
cussed in  [Grad  1964a,  1964b,  1966]. 
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4 .   Large  Mean  Free  Path- -The  Qbiquitous  Logarithm 

In  the  limit  of  large  mean  free  path  we  write  e  =  K 
(Knudsen  number).   There  are  two  distinct  reasons  why 
K  — >  00  does  not  approach  uniformly  to  the  free  flow  limit. 
One  reason  is  that  ^(^)  — ^0  for  small  i   so  there  is  always 
a  class  of  slow  particles  with  small  free  path  no  matter 
how  large  K  may  be.   A  second  reason  is  that  in  an  exterior 
or  otherwise  unbounded  domain  there  will  always  be  distant 
collisions . 

There  is  a  simple  mathematical  lemma  which,  although 
mathematically  trivial,  is  responsible  for  many  important 
phenomena  which  are  physically  profound  and  subtle,  includ- 
ing both  of  the  above  mentioned  sources  of  nonunif ormity 
at  large  K.   The  divergent  integral  jdr/r  can  be  given  a 
finite  value  by  introducing  a  convergence  factor 


Ke) 


0(er)  ^  (4.1) 


(the  divergence  near  r  =  0,  which  we  dismiss  with  the 
finite  limit  r  =  1,  can  be  studied  in  exactly  the  same 
way).   We  assume  that  0{O)    =   1   and  52f(s)  is  smooth  near 
s  =  0;  also  that  (Zf  is  integrable  at  infinity.   Writing 
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p(3)f  =J(^-1)?^| 


1 

ds_ 
s 
e 


(/>   ^        (4.2) 
"^  s 


we  see  that  the  first  and  third  integrals  are  finite 
as  6  — >  0.   Consequently, 

1(e)  =  log  1/e  +  0(1)  .  (4.3) 

We  observe  that  the  dominant  term  log  l/e  is  completely 
independent  of  the  nature  of  the  convergence  factor,  depend- 
ing only  on  the  length  scale  e;  e.g.,  its  value  would  be  the 
same  if  e  were  replaced  by  2e.   Also  the  dominant  contribu- 
tion to  the  Integral  comes  from  values  of  r  which  are  large 
compared  to  the  lower  limit,  r  =  1,  but  are  small  compared 
to  r  =  I/e  [where  0"  deviates  from  the  value  0(0)  =  1]. 
The  error  term,  0(1)  in  (4.3),  does  depend  on  the  entire 
function  0( s ) .   In  other  words,  the  leading  term  of  an 
expansion  of  1(e)  is  local  and  extremely  insensitive, 
whereas  higher  order  terms  are  global  and  more  complicated. 
The  special  choice 


0(s)=l,   0<s<a 

(4.4) 
0(s)  =  0  ,   s  >  a 
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'r 


which  amounts  to  cutting  off  the  simple  integral   dr/i 
at  r  =  a/e,  also  gives  a  leading  term  which  is  independent 
of  the  cutoff  parameter  a. 

An  integral  of  the  form  (4.1)  arises  commonly  when  a 
problem  presents  two  length  scales,  as  in  a  boundary  layer 
analysis.   More  generally. 


CO 

J(e)  =  [  0{ev)   ^  (4.5) 


can  be  expanded  as 

(4.6) 
There  is  a  polynomial  in  e  ending  with  the  appearance  of 
a  logarithm.   All  the  listed  coefficients  are  local 
[depending  on  derivatives  of  j25  at  r  =  0];  higher  order  terms 
In  the  expansion  depend  globally  on  j6. 
Similar  expansions  can  be  made  for 

00 

K(e)  =  J  0(er)V'(r)dr  (4.7) 

1 

In  terms  of  as  many  moments  of  ^(r)  as  are  finite  and  an 
equivalent  number  of  derivatives  of  0   at  the  origin.   The 
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boundary  layer  interpretation  of  this  result  states  that 
the  dominant  part  of  the  expansion  of  K{e)  (up  to  and 
including  the  logarithm)  can  be  evaluated  locally,  within 
the  boundary  layer,  without  the  necessity  of  matching  to 
the  exterior  region. 

A  common  tool  in  investigating  solutions  of  the 
Boltzmann  equation  is  to  write  it  as  an  integral  equation. 
In  the  integrand  will  appear  a  factor 

exp[-v(^)x/^^]  =  exp[-x/A(U]  •  (^.8) 

Near  a  wall,  x  ^  0,  we  can  expand  the  exponent  in  powers 
of  X  except  for  small  ^-,  where  A  «,  0.   Integrals  with 
respect  to  ^  of  such  an  expression  will  yield  a  leading 

term  x  log  l/x  (or  higher  order  terms  x  +  x   +  ... 

r     / 
+  x  log  l/x  for  moments  which  have  powers  of  ^,  as  a 

factor) . 

For  large  mean  free  path,  the  relevant  factor  appears 

as 

exp[-x/0(^)]  (4.9) 

With  X  fixed  and  nonzero,  this  "convergence  factor"  is 
regular  in  1/k  except  where  A  ^  0.   Again  we  obtain 
(l/K)log  K  as  the  leading  term  for  an  integral  over  this 
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expression.   But  this  Is  correct  only  In  one  dimension. 

At  a  point  definitely  removed  from  a  wall,  the  dimension 

of  the  manifold  on  which  the  appropriately  defined  A(4) 

vanishes  depends  on  the  dimensionality  of  the  problem. 

In  one  dimension  (e.g.,  two  parallel  flat  plates),  A  =  0 

when  ^,  =  0  and  (^p,^)  are  arbitrary.   In  two  dimensions 

(e.g.,  two  concentric  cylinders),  A  =  0  when  ^  =   ^a  =   0 

r    y 

and  ^   is  arbitrary.   In  three  dimensions  (e.g.,  two 
concentric  spheres,  or  any  bounded  interior  domain), 
A  =  0  only  when  the  vector  |  =  0.   Therefore,  in  two 
dimensions  there  will  be  a  term  1/K  followed  by  (log  K)/K  , 
and  in  three  dimensions  two  terms,  1/k  +  1/k  ,  followed  by 
(log  K)/K  .   Furthermore,  the  evaluation  of  these  coeffi- 
cients is  local,  depending  essentially  on  the  free  flow 
distribution,  whereas  any  higher  order  terms  pose  diffi- 
culties of  matching  in  velocity  space. 

In  an  external  domain  we  have,  in  addition  to  the  slow 
particles  treated  above,  the  effect  of  distant  particles. 
The  same  integral  (4.9)  appears  as  a  convergence  factor, 
now  integrated  with  respect  to  x  from  Infinity  Instead  of 
from  a  wall.   In  two  dimensions  (e.g.,  flow  past  a  cylinder), 
the  function  which  is  multiplied  by  this  exponential 
convergence  factor  decays  as  l/x  (this  represents  the  free 
flow  disturbance  which  is  proportional  to  the  angle 
subtended  by  the  object  at  a  distance  x).   Thus  the  leading 
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term  is  log  K.   An  additional  factor  1/K  Is  explicit  in 
the  Integral  which  is  evaluated;  thus  the  modification 

of  the  free  flow  drag  is  (log  K)/K.   Similarly,  in  three 

o 
dimensions  where  the  influence  function  decays  as  l/x  , 

we  obtain  1/K  +  (log  K)/K  as  the  leading  terms  [after 

free  flow  which  is  scaled  to  be  0(1)]. 

In  an  exterior  domain  where  both  the  slow  and  distant 
particle  nonunif ormities  are  present  (arising  from  inte- 
gration over  ^  and  x  respectively),  it  is  the  distant 
particle  effect  which  dominates,  always  entering  at 
lower  order.   Alternatively,  the  logarithm  enters  at  a 
lower  order  in  a  two-dimensional  exterior  problem,  say 
drag  past  a  cylinder,  than  in  a  two-dimensional  Interior 
problem,  such  as  drag  between  two  concentric  cylinders. 

We  have  been  intentionally  vague  in  not  stating 
above  exactly  what  integral  is  being  evaluated.   The 
descriptions  given  are  basic  in  the  sense  that  there 
exists  some  Interesting  macroscopic  quantity  with  the 
logarithm  appearing  at  the  given  order.   Other  observ- 
ables,  involving  integrals  over  functions  which  vanish 
at  small  velocity  or  Involving  some  special  symmetry, 
may  exhibit  a  logarithm  only  later  on  in  the  expansion. 
For  example.  In  linearized  one-dimensional  heat  flow, 
although  the  temperature  gradient  at  the  midplane  is 
(log  K)/K   [Willis,  1958],  the  ratio  of  total  heat  flow 
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to  the  difference  In  wall  temperature  (effective  conduc- 
tivity) is  1+I/K  (because  the  heat  flow  is  a  moment  of  f 
with  a  factor  ^-|  ) .   But  in  a  nonlinear  problem  with  a 
finite  temperature  difference,  the  effective  conductivity 
(normalized  at  fixed  total  mass  in  the  system)  will  behave 
like  1  + (log  K)/K  because  of  a  loss  of  symmetry  about  the 
midplane. 

The  integrals  which  arise  in  this  subject  will,  in 
general.  Involve  combinations  of  boundary  layers  x  •,  0, 
^  «,  0  and  nonunif ormities  x  ^  <»,  K  «.  <».   Although  the  basic 
concepts  and  mathematical  analysis  of  these  linked  non- 
uniformities  are  elementary,  the  bookkeeping  can  be  quite 
complicated,  involving  the  number  of  dimensions,  the  degree 
of  symmetry,  and  the  specific  observable.   A  complete  list 
has  not  been  written  down,  although  many  specific  examples 
have  been  studied,  e.g.  [Willis,  1958;  Cercignani,  19^2; 
Su,  1964;  Sone,  1964;  Pao,  I965;  Grad,  1966;  Yu,  1967]. 

It  is  interesting  to  note  tha c  che  same  mathematical 
lemma,  and  approximately  the  same  physical  situation, 
occurs  in  the  derivation  of  a  kinetic  equation.   To  lowest 
order  in  the  appropriate  expansion  parameter  through  which 
the  Boltzmann  equation  arises,  there  is  no  formal  difficulty. 
To  higher  order,  including  dense  gas  corrections,  one  finds 
a  nonunif ormity  arising  from  slow  particles  and  another  from 
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distant  particles  [Grad,  1958;  Prieman  and  Goldman,  I966, 
and  many  others ] . 

Essentially  the  same  mathematical  situation  also 
arises  in  evaluating  transport  coefficients  in  a  plasma. 
A  formal  divergence,  arising  from  the  long  range  Coulomb 
potential,  is  eliminated  on  introducing  the  concept  of  a 
shielded  potential  [Spitzer,  I956] .   Without  performing 
any  of  the  lengthy  computations,  it  is  immediately  evident 
that  the  dominant  term  of  the  electrical  conductivity 
(i.e.,  the  Coulomb  logarithm)  must  arise  with  identical 
numerical  coefficients  whether  we  employ  the  Fokker-Planck 
equation  or  the  Boltzmann  equation,  and  whether  we  use  a 
spherical  Debye  shielding  cloud,  or  a  more  sophisticated 
distorted  shielding  cloud,  or  whether  we  simply  cut  off 
the  divergence  at  the  Debye  length  [Grad,  1962].   The  only 
fact  that  we  must  know  is  the  scale  length  of  the  shielding, 
without  regard  to  the  details  of  the  shielding  process. 
But  to  carry  the  computation  beyond  the  dominant  term  Is 
formidable  and  requires  a  correct  model  for  the  shielding. 
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5 •   Small  Mean  Free  Path- -Introduction 

The  limit  of  small  mean  free  path  probably  represents 
the  most  extreme  of  all  nonunif ornnities  to  be  found  within 
the  scope  of  the  Boltzmann  equation,  but  it  is  also  the 
best  understood  and  documented.   In  this  limit,  the  gas 
description  is  reduced  from  a  distribution  function 
f(^,x,t)  to  the  macroscopic  fluid  coordinates  p(x,t), 
T(x,t),  u(x,t). 

The  Boltzmann  equation  can  be  written  in  skeleton  fonn 
as 


||+l-||=Q(f,f)  (5.1) 


where  Q  is  a  quadratic  operator  on  f  which  vanishes  when 
applied  to  a  local  Maxwellian 


^M  ^  S72   ^^Pf-  4^  ]  •  (5.2) 


By  local  we  mean  that  p,  T,  and  u  may  be  arbitrary  functions 
of  X  and  t.   Under  appropriate  scaling,  the  Boltzmann 
equation  becomes 


-  26  - 


If +4  •  ||  =  i«(f,f)  (5.J) 


or  In  a  linearized  version 


^+  l-||  +  il.(f)  =0  (5.*) 


The  nonlinear  form  can  be  obtained  by  the  substitution 
f  — >  f /e  In  {5«1)  and  the  linear  version  by  the  substi- 
tution f  — >  —  f j^  +  f  where  f„  Is  here  an  absolute  Max- 
welllan  with  constant  parameters  p,  T,  u.   The  signifi- 
cance of  the  substitution  f  — )  f/e  is  that  for  small  e 
the  density  Is  large  and  therefore  the  mean  free  path  is 
small.   In  the  linearization,  e  represents  the  mean  free 
path  of  the  unperturbed  constant  state,  f,y,/S' 

It  is  evident  that  the  limit  e  — >  0  is  singular  since 
all  derivatives  are  lost.   For  small  e  we  should  expect  to 
have  Q(f,f )  „  0  and  f  «  f   except  in  regions  where  ^/^t   or 
h/hx   become  unbounded.   To  be  more  definite,  take  a  fixed 
domain  with  appropriate  fixed  boundary  and  initial  condi- 
tions on  f,  independent  of  the  parameter  e.   The  qualita- 
tive behavior  for  small  e  of  the  solution  f(^,x,t;  e)  is 
described  by  a  conjecture  [Grad,  1958]  which  states  that  f 
approximates  f„,  and  the  fluid  parameters  in  f„  obey  fluid 
equations,  except  where  this  is  evidently  impossible. 
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Specifically,  there  are  three  locations  where  f  cannot 
approximate  f„  closely: 

1)  In  an  Initial  layer,  since  the  initial  value 
f(^,x,0)  is  arbitrary; 

2)  in  a  boundary  layer  where  the  "half"  distribu- 
tion of  molecules  leaving  the  wall  is,  in  principle, 
arbitrary; 

3)  in  a  shock  layer  where  one  cannot  conclude  that 
Q(f,f)  «  0  since  S/^t  and  b/hx   are  of  order  l/e  (the  shock 
thickness  scales  as  e). 

The  conjecture  states  that,  except  in  these  three 
layers  of  nonunlf ormity,  the  Hilbert  or  Chapman-Enskog 
formal  expansions  in  e  are  asymptotic  to  the  actual  solution 
f(^,x,t;  e).   This  conjecture  has  been  promoted  to  a  proved 
theorem  in  linear  and  weakly  nonlinear  formulations  of  the 
initial  value  problem  [Grad,  1963a, b;  I965]. 

Accepting  the  conjecture  more  generally,  the  question 
immediately  arises  of  matching  the  asymptotic  expansions 
across  the  three  types  of  boundary  layers  within  which  they 
are  not  valid.   The  fluid  equations  (which  can  take  several 
forms)  yield  solutions  only  when  they  are  completed  with 
initial  data,  boundary  data,  and  possibly  shock  jump  condi- 
tions.  The  relation  between  the  given  Boltzmann  boundary 
data  (half  distribution)  and  the  proper  asymptotic  fluid 
boundary  data  requires  analysis  of  a  boundary  layer,  and 
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similarly  for  the  Initial  and  shock  layers.   A  crude  anal- 
ysis of  the  boundary  layer  was  given  by  Maxwell  and  led 
to  the  concept  of  boundary  slip  (of  tangential  velocity 
and  temperature).   In  the  same  terminology  we  can  talk  of 
Initial  slip  and  shock  slip.   The  initial  slip  represents 
the  amount  by  which  the  true  initial  values  of  the  fluid 
state  (given  as  moments  of  the  assigned  initial  distribu- 
tion function)  must  be  altered  so  that,  when  used  as 
initial  values  In  fluid  equations,  they  give  rise  to  the 
correct  future  fluid  state  as  it  would  be  found  by  solving 
the  Boltzmann  equation.   Similarly,  the  shock  slip  repre- 
sents the  correction  to  the  classical  Hugoniot  condition 
which  is  required  for  a  fluid  calculation  to  yield  the 
same  result  as  solution  of  the  Boltzmann  equation.   A 
related  concept  of  shock  slip  has  also  been  found  within 
fluid  dynamics  proper  as  a  modification  of  the  Invlscld 
Hugoniot  conditions  to  take  into  account  the  interaction 
between  a  shock  and  a  boundary  layer  [Probstein  and  Pan, 
1963].   Although  viscous  fluid  equations  can  presiimably 
give  a  smooth  shock  transition,  to  approximate  the  kinetic 
solution  correctly,  they  must  be  allowed  to  have  disconti- 
nuities . 

The  small  mean  free  path  limit  has  been  split  into 
two  problems.   The  first,  in  the  interior  (away  from 
Initial,  boundary,  and  shock  layers),  is  solved  by  a 
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Hllbert  or  Chapman -Ens kog  expansion  (except,  as  we  shall 
see  later,  when  they  require  modification).   The  second 
consists  in  the  analysis  of  three  kinetic  boundary  layers 
In  order  to  complete  the  fluid  formulation  with  initial, 
boundary,  and  jump  conditions.  Note  that  the  type  of 
boundary  data  that  is  required  for  invlscld  fluid  equa- 
tions is  quite  different  from  what  is  required  for  viscous 
equations.   The  kinetic  boundary  layer  analysis  will  there- 
fore differ  depending  on  what  fluid  model  is  selected. 
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6.   The  Initial  Value  Problem 

The  solution  to  the  problem  of  Initial  slip  Is  essen- 
tially complete  [Grad,  I963;  1965]*   For  the  nonlinear 
problem  the  solution  is  formal;  an  algorithm  is  given  for 
the  computation  of  the  asymptotic  initial  fluid  state  in 
terms  of  the  given  initial  distribution  function.   For 
the  linear  and  weakly  nonlinear  problems,  the  solution 
f(^,x,t;  e)  has  been  shown  to  exist,  and  the  formal  expan- 
sions in  £  have  been  proved  to  be  asymptotic  to  this  solu- 
tion. 

To  be  more  precise,  let  f(^,x,t;  e)  be  the  solution 
of  the  Boltzmann  equation  which  takes  the  initial  value 
f(^,x,0).   Prom  f(^,x,0),  we  extract  its  moments  p(x,0), 
u(x,0),  T(x,0)  and  solve  the  inviscid  Euler  equations  with 
these  initial  values.   Let  f„(|,x,t)  be  the  local  Maxwelllan 
formed  from  these  Euler  solutions  p(x,t),  u(x,t),  T(x,t). 
Then 


lim  f(^,x,t;  e)  =  f„(^,x,t)  (6.1) 

e  -»  0  ^ 


for  any  fixed  t  >  0  and  in  a  variety  of  norms  in  (x,4) 
(Including  Lp,  maximum,  and  others).   In  the  limit  e  — >0 
the  Boltzmann  equation  disappears;  but  its  solutions  have 
a  well  defined  limit.   The  moments  p,  T,  u  are  continuous  at 
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t  =  0  uniformly  In  e^  but  f  Is  not;  nor  are  ^u/Bt  and  ST/St 

The  Navier-Stokes  equations  must  be  supplied  with 
corrected  initial  data,  u  (x,0;  e)  =  u  (x,0)  +  eu-,(x,0), 
etc   Here  Up^(x,0)  is  the  actual  initial  velocity  [a 
moment  of  f(|,x,0)],  and  the  correction  u, (x,0)  is  ob- 
tained from  f(^,x,0)  by  a  definite  analytic  procedure. 
The  solution  p  (x,t;  e),  T  (x,t;  e) ,    u  (x,tj  e)  of  the 
Navier-Stokes  equations  (in  which  e  multiplies  the  vis- 
cosity and  heat  conductivity)  is  used  to  construct  a 
Chapman-Ens kog  distribution  function  f^  =  f „  +  ef,  where 
f   is  locally  Maxwellian  with  p  ,  Tj^,  u  as  parameters, 
and  f-,  depends  on  gradients  of  u  and  Tj^.   It  is  then 
shown  (in  various  norms)  that 


lim  ^  If  -  f  1  =  0  ,  (6.2) 

e  ->  0  ^       ^ 


this  time  uniformly  in  t  >  t_ . 

The  quantities  p-,,  T-^,    u,,  represent  the  initial  slip, 
i.e.,  the  discrepancy  between  the  initial  fluid  state  and 
the  initial  value  which  is  required  for  the  Navier-Stokes 
solution  to  correctly  approximate  the  Boltzmann  solution. 
The  first  order  initial  slip  vanishes  if  f(|,x,0)  is  a 
local  Maxwellian.   The  initial  slip  vanishes  to  higher 
order  if  the  initial  distribution  is  a  Chapman -Ens kog 
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function;  but  this  requires  the  initial  function  f(^,x,0;  e) 
to  depend  explicitly  on  e. 

Proofs  of  these  statements  make  use  of  strong  exis- 
tence theorems  with  estimates  of  f  and  of  space  derivatives 
of  f  that  are  uniform  in  the  singular  parameter  e. 

The  statement  that  a  Boltzmann  solution  is  approxi- 
mated by  fluid  equations  outside  the  regions  of  nonuni- 
formity  is  somewhat  flexible.   The  Hilbert  theory  (linear 
or  nonlinear)  gives  a  formal  expansion  of  f  in  powers  of  e 
governed  by  the  invlscid  Euler  equations  to  lowest  order 
and  by  successive  inhomogeneous  Euler  equations  to  higher 
order.   More  precisely,  the  higher  order  equations  are  the 
inhomogeneous  variational  equations  of  the  Euler  system 
(linear  with  variable  coefficients,  having  the  same  char- 
acteristics as  the  Euler  equations).   The  Chapman- Ens kog 
theory  also  gives  the  Euler  equations  to  lowest  order; 
but  this  is  followed  by  the  Navier-Stokes  equations,  the 
Burnett  equations,  and  a  sequence  of  partial  differential 
equations  in  p,  T,  and  u  of  successively  higher  order  in 
space  derivatives.   In  other  words,  the  Chapman-Enskog 
theory  gives  a  formal  expansion  of  the  macroscopic 
equations  which  govern  the  fluid's  evolution,  whereas  the 
Hilbert  theory  gives  a  formal  expansion  of  the  solution. 
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Both  expansions,  Hllbert  and  Chapman-Enskog,  are 
asymptotic.   But  the  ranges  of  validity  are  not  the  same. 
Each  expansion,  when  truncated,  satisfies  the  Boltzmann 
equation  with  a  remainder. 


-^  +    i    '    ^  +  -  Uf    )    =    s'r   (|,x,t;  e)  (6.3) 

dt        ox    en       n  ^ 


The  Hllbert  remainder,  which  is  obtained  by  solving  inhomo- 
geneous  Euler  equations,  grows  algebraically  as  t  .   The 
Navler-Stokes  remainder  decays  as  t  — >  oo.   This  is  why  the 
Hllbert  expansion  is  asymptotic  for  only  limited  time, 
tp^  <  t  <  t-,  ,  whereas  the  Navler-Stokes  result  is  uniformly 
valid  for  t  >  t_.   The  higher  order  Chapman-Enskog  approxi- 
mations, Burnett,  etc.,  are  only  formal  at  the  present  time 
in  the  absence  of  any  estimates  of  the  boundedness,  growth, 
or  decay  of  their  solutions. 

The  proof  that  a  given  formal  expansion  is  asymptotic 
to  a  solution  of  the  Boltzmann  equation  is  quite  oblivious 
to  the  recipe  which  gave  rise  to  the  expansion.   All  that 
needs  to  be  verified  is  that  the  expansion  satisfies  the 
Boltzmann  equation  term  by  term  with  a  remainder 
e'^R  (|,x,t;  e)  where  R  is  suitably  bounded.   ¥e  may  con- 
struct such  a  series  by  ad  hoc  means  as  well  as  by  expan- 
sion.  For  example,  from  the  Hllbert  or  the  Chapman-Enskog 
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formulas  we  can  easily  construct  a  unique  sequence  starting 
with  Euler,  Navler-Stokes,  and  continuing  with  Inhomogeneous 
Navler-Stokes  equations.   It  Is  trivial  to  verify  that  this 
sequence  is  asymptotically  equivalent  to  either  the  Hilbert 
or  Chapman-Ens kog  series.   This  formal  series  satisfies 
the  Boltzmann  equation  with  an  error  term  R  which  decays 
for  large  t  at  any  order  n  >  1.   This  series  is  therefore 
superior  to  the  Hilbert  series  in  range  of  validity,  and  it 
is  superior  to  the  Chapman-Enskog  series  since  the  theory 
of  the  associated  partial  differential  equations  is  known. 
Although  this  procedure  is  arbitrary  for  the  Initial  value 
problem,  we  shall  find  that  it  arises  naturally  in  slow 
steady  flow  (Sec.  8). 

The  solution  of  the  initial  slip  problem  can  be  given 
a  simple  abstract  formulation  in  the  case  of  the  linear 
Boltzmann  equation  (5-^)-   In  terms  of  a  suitable  normaliza- 
tion, the  local  Maxwellian  takes  the  form 


1  1   ? 

00  =  ^7^  exp{-  -p  ^  )  (6.4) 

(2Tr)^/^        ^ 


and  the  nullspace  of  the  collision  operator  L, 


L(^^)  =  0  ,    a  =  0,1,2,3,4  (6.5) 
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Is   given  by 


^0  =  03^/2 

^1  ^   ^X^^  i  =  1'2,3  (6.6) 


^4  -   (^^  -  3)(co/6)l/2 


The  i}/     are   orthonormal    In   the   L„  norm. 


(f 


,g)   =J   f  g  d?  (6.7) 


The  (linearized)  fluid  state  Is  given  by  the  projection 
of  f  Into  the  nullspace  of  L, 

f  =y   9  ii  (6.8) 


where 


(^o'f)  =  Po  =  P 

(^^,f)  =  p^  =  u^    1  =  1,2,3  (6.9) 

(^i^.f)  =  P4  =  (5/2)^/^T 


We  Introduce 


f  +  f  (6.10) 
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where  f  is  the  projection  into  the  orthogonal  complement 
to  the  nullspace. 

The  basis  of  the  Hilbert  theory  is  the  hypothesis 
that  f(^,x,t;  e)  be  expressible  as  a  formal  power  series 
in  e.   The  essential  conclusion  is  that  to  every  given 
fluid  state,  p  (x),  is  assigned  a  unique  distribution  func- 
tion [Hilbert,  1912] 


fjj(e,x;  e)  =  f°  +  efjj  +  ...  (6.11) 


An  equivalent  statement  is  that  to  a  given  projection 
f(^,x),  there  is  assigned  a  unique  distribution  f„;  in 

n 

particular,  the  leading  term  f„  is  exactly  the  given 

n 

fluid  state  f,  and  all  succeeding  terms  are  in  the 
orthogonal  complement.   The  manifold  of  special  func- 
tions f„  is  formally  invariant  under  solution  of  the 
n 

Boltzmann  equation.   In  other  words,  a  special  function 
f„  taken  as  an  initial  value  for  the  Boltzmann  equation 

ri 

will  give  rise  to  a  solution  which  remains  within  the 
space  of  functions  f..-   A  complementary  result  has  also 

n 

been  proved  [Grad,  1965a].   The  hypothesis  that  f(^,x,T;  e), 
(where  x  =  t/e)   is  a  power   series  in  e  and  that  f 
decays  to  zero  implies  that  to  a  given  projection  F(^,x) 
is  assigned  a  unique  distribution  function 
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f^  =  f°  +  ef^  +  ...  (6.12) 


where  f^  Is  Just  the  given  projection  f ,  and  all  succeeding 
terms  lie  in  the  nullspace  (i.e.,  they  are  strictly  fluid 
components).   This  procedure  also  defines  an  asymptotically 
Invariant  space  of  solutions  of  the  Boltzmann  equation. 

These  two  results  can  be  combined  to  obtain  the  unique 
resolution 


f(?,x)  =  r^(i,x;    e)    +  fg(^,x;  e)  (6.13) 


of  an  arbitrary  distribution  f,  which  we  choose  to  be 
independent  of  e.   To  simplify  the  notation  we  write  h 
for  f„  and  g  for  f_.   The  orthogonal  decomposition 

rl  (j 

f  =  F  +  f  of  the  given  function  f(i,x)    is  used  to  obtain 
the  leading  terms 


hg  =  f  ,   gg  =  f  (6.ii^) 


which  generate   the   two  series 


hQ(  ^,x;    e)    =   h^   +   eh^   +    .  . 


gQ(C,x;    e)    =  gQ   +   egQ   +    .. 


(6.15) 
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The  difference 

f  -  hQ  -  gQ  =  ef^(e,x;  e)  (6.16) 

>\ 

is  similarly  decomposed,  ^1=^1"'"  ^1  ^°  generate  the 

two  series  h,  and  g, ,  etc.   In  this  way  we  obtain  (6.13) 
as  the  series 


f  =  f  -  +  eh,  +  .  .  . 
H    0     1 


^G  =  So  +  Bgi  +  .  . 


(6.17) 


Let  f(C,x,t;  e)  be  the  solution  of  the  Boltzmann 
equation  with  f(t,x)  as  Initial  value.   Taking  as  an 
initial  value  for  the  Boltzmann  equation  the  Hllbert 
component  f„  of  the  resolution  (6. 13)  of  f(^,x)  yields  a 

n 

solution  which  Is  asymptotically  equal  to  f(^,x,t;  e)  after 
a  time  0(e).   The  complementary  solution,  taking  f^  as 
initial  value,  describes  a  transient.   Although  the  time 
evolution  of  f..  Is  governed  by  fluid  equations,  the  entire 

rl 

initial  function  f(^,x)  is  reflected  in  the  fluid  initial 
data,  f„(^,x;  e). 

If  the  collision  operator  L  has  discrete  eigenvalues 
other  than  zero,  then  the  "transient"  component  f-  can 
be  further  decomposed  into  fluidlike  components  which 
satisfy  partial  differential  equations  (similar  to  Euler, 
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Navler-Stokes,  etc.)?  but  which  decay  exponentially  in 

time. 

If  f{i,x)   has  n  derivatives  with  respect  to  x,  the 

separation  into  f„  and  f~  can  be  carried  out  only  to 

n       u 

order  e"^"  .   The  problem  of  discontinuous  or  more  general 

* 
non-smooth  initial  data  remains  an  open  question.   For 

Maxwell  molecules,  the  expansion  f^  is  explicitly  conver- 
gent, as  a  polynomial  in  e,  when  f(4,x)  is  a  polynomial 
in  X  [Grad,  1963a].  One  can  expect  convergence  of  the 
series  for  f„  if  f  is  sufficiently  smooth  (analyticlty  is 

n 

not  enough).   A  conjecture  made  in  [Grad,  1963a],  based  on 
Maxwellian  molecules,  is  that  convergence  requires  f  to 
be  an  entire  function  of  x  of  order  one-half.   In  the 
special  case  of  hard  spheres,  a  single  Fourier  coefficient, 
f  =  0(^)e    ,  has  been  shown  to  give  rise  to  a  convergent 
series  for  f„  when  k  is  sufficiently  small  [McLennan,  I965] 
It  is  not  clear  whether  an  entire  function  of  order  one 
(like  e    )  is  sufficient  for  convergence  in  general,  or 
whether  there  is  an  essential  difference  between  hard 
spheres  and  Maxwellian  molecules. 


For  a  model  of  the  Boltzmann  equation,  Y.  Sone  [Phys. 
Fluids,  to  appear]  has  obtained  the  asymptotic  fluid 
behavior  for  a  special  class  of  discontinuous  initial  data. 
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7'   The  Boundary  Value  Problem- -General  Discussion 

The  boundary  layer  presents  a  more  difficult  problem 
than  the  initial  layer  for  at  least  two  reasons.   First  of 
all,  the  basic  proof  of  existence  is  more  difficult.   There 
have  been  great  strides  recently  in  existence  theory  for 
the  boundary  value  problem,   but  these  results  are  not  yet 
adequate  to  the  demands  of  the  small  mean  free  path  limit. 
Neither  the  required  uniform  bounds  as  e  — )  0  nor  estimates 
of  derivatives  are  yet  available.   This  is  not  simply  a 
technicality,  since  space  derivatives  do  become  infinite 
near  a  boundary. 

The  second  difficulty  is  that  the  limit  e  — >  0  is 
ambiguous  even  in  fluid  dynamics.   We  recall  that  the 
viscosity  and  heat  conductivity  coefficients  are  propor- 
tional to  e.   In  a  Stokes  linearization  about  a  free  stream 
U  =  0,  the  flow  is  essentially  independent  of  value  of  the 
viscosity  coefficient;  cf.  slow  flow  about  a  sphere  or 
Poisseullle  flow  with  a  parabolic  velocity  profile.   On 
the  other  hand,  in  an  Oseen  linearization  about  a  finite 


C.  Cercignani  (private  communication)  in  a  general 
bounded  domain  for  the  linear  equation  and  Y.  P.  Pao 
(private  communication)  for  the  weakly  nonlinear  equation 
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velocity  U,  the  limit  \i   — >  0  is  singular  and  gives  an 

inviscid  limiting  flow  except  within  a  boundary  layer 

1/2 
and  wake.   This  boundary  layer,  of  thickness  0( e  ^  ), 

is  distinct  from  the  kinetic  boundary  layer  which  is  0(e). 

To  obtain  a  definite  result  in  the 'kinetic  problem, 
we  must  at  least  distinguish  between  the  various  macro- 
scopic cases.   One  distinction  is  the  value  of  the  Reynolds 
number.   In  the  Stokes  linearization  we  have  set  Re  =  0 
by  setting  the  free  stream  velocity  equal  to  zero  before 
[i   — >0.   In  the  Oseen  linearization  Re  is  finite,  but  it 
becomes  Infinite  in  the  subsequent  limit  \x   — >0.   We  shall 
find  it  convenient  to  Interpolate  between  these  cases  by 
taking  an  Oseen  linearization  after  which  we  allow  U  to 
approach  zero  with  e.   In  this  limit  Re  remains  finite. 

The  interior  of  the  fluid  will  be  described  by  a 
formal  expansion  which  is  equivalent  to  Hilbert's  in  the 
initial  value  problem.   By  supplementing  this  expansion 
with  conjectured  existence  theorems,  we  shall  also  be 
able  to  imitate  the  rigorous  analysis  of  the  initial 
layer  by  a  plausible  analysis  of  the  boundary  layer.-   This 
will  give  the  asymptotic  boundary  conditions  which  are 
needed  to  complete  the  formal  fluid  expansion  in  the 
interior. 

In  the  case  of  the  Stokes  linearization  and  the  more 
general  finite  Re  linearization  just  described,  we  find 
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that  the  formal  Hilbert  expansion  automatically  gives 
Navier-Stokes  followed  by  successive  Inhomogeneous 
Navler-Stokes  equations,  rather  than  Euler  and  Inhomo- 
geneous Euler  equations  as  in  the  initial  value  problem. 
The  reason  is  that  the  steady  flow  Euler  equations  are 
too  degenerate.   They  are  not  a  determined  system  of 
differential  equations  which  yield  solutions  merely  by 
supplementation  with  boundary  conditions.   In  common  with 
most  asymptotic  expansions,  a  given  term  in  the  standard 
Hilbert  expansion  is  not  uniquely  determined  until  we  turn 
to  the  next  higher  order  from  which  we  extract  compati- 
bility conditions.   In  the  case  of  steady  flow,  the 
additional  degeneracy  requires  us  to  extract  information 
from  two  succeeding  levels  to  complete  the  determination 
of  a  given  order.   Thus  the  expansion  of  a  solution  of 
the  Boltzmann  equation  in  powers  of  e  is  essentially 
different  in  steady  flow  and  time-dependent  flow. 

In  the  case  Re  — )  <»,  a  restriction  to  the  viscous 
boundary  layer  rather  than  the  entire  Interior  allows 
a  similar  result  (inhomogeneous  Navler-Stokes  equations) 
to  be  obtained  [Trilling,  1964]. 

For  the  general  case,  nonlinear  and  nonsteady,  we  can 
find  an  asymptotic  fluid  description  by  an  ad  hoc  proce- 
dure similar  to  the  one  described  above  (Sec  6)  in  which 


^3 


we  proposed  Inhomogeneous  Navler-Stokes  equations  for 
the  initial  value  problem.   The  Important  point  with 
regard  to  the  boundary  layer  Is  that  the  boundary  condi- 
tions which  it  must  supply  are,  to  any  order  in  e,    Just 
what  Is  appropriate  to  a  Navler-Stokes  system. 
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8.   steady  Hilbert  Expansion  for  Finite  Reynolds  Number 
¥e  take  the  linear  Boltzmann  equation 


11+  ?-||  +  iL(f)  =0  (8.1) 


in  the  normalized  and  dimensionless  form  described  in 
equations  (6.4)  -  (6.9).   A  Stokes  linearization  in 
steady  flow  is  simply 


br    .    1 
e 


i'  ^  +  T   L(f)  =  0  (8.2) 


and  an  Oseen  linearization  about  the  free  stream  velocity 
eU  yields 


^"•i^  «-||^iL(f)  =0  (8.3) 


(the  molecular  velocity  ^  is  relative  to  the  stream  at 
infinity,  as  is  appropriate  in  taking  moments  of  f ) .   We 
note  the  linearized  equation  of  state 

p  =  p  +  T  (8.4) 

and  the  conservation  equations  [appropriate  to  (8.1)] 
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If  +  dlv  u  =  0 


-^-^-^-°  (8-5) 


■^  (|-  T)  +  dlv  u  +  div  q  =  0 


The  dlmenslonless  stress  devlator  and  heat  flow  vector 
are  the  projections 


P,j  =  (^,^j  -J^%J'  ^) 


^i  =  (|  f^^  -  5]?!^  f) 


(8.6) 


The  Euler  equations  are  obtained  by  setting  p.  .  =  q.  =  0 
In  (8.5) J  and  the  Navler-Stokes  equations  result  from 


fSu .    Su  .   p  ") 

3— i  +  -^— i  -  X  5.  .  dlv  u  y 
dx  .    dXj_   5   ij       J 

(8.7) 

St 


^^-  =  -^  ^ 


'1 


The  full  Chapman- Ens kog  expansion  takes  the  form  [Grad,  1963a] 
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|£-  +  dlv  u  -  0 


3u  ,  „_   ''\~  _2n+l  .n^ 

n     "n 


^  +  Vp  =  ^2_  e     A  [m.„Au  +  a^  V(dlv  u)] 


0 


00 


+  2_  e^"^  ^'^[Pn  VT  +  7^  Vp] 


(8.8) 


00 


^  (|  T)  +  dlv  u  =  i:  e^^^^   A^^l  [>^T  +  5^p] 


+  2_  e^"  ^n  A'^ldlv  u) 
1       " 


where  the  coefficients  |j,,X,a,p,7,5,C   depend 

n   n   n   n   'n   n   n 

on  the  molecular  model.   In  particular,  |i  =  e|i  Is  the 
(dlmenslonless )  viscosity  coefficient  and  A  =  eA   is  the 
heat  conductivity  (a  =  11/3  and  5_  =  0).   The  steady 
flow  Stokes  linearization  is  obtained  by  setting  h/ht  =   0, 
and  the  Oseen  linearization  is  obtained  by  replacing  B/St 
by  U  •  V. 

In  a  steady  Stokes  flow  the  Euler  equations  are 


dlv  u  =  0 

Vp  =  0  (8.9) 
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and  the  Navler-Stokes  equations  are 

dlv  u  =  0 

Vp  +  M,  curl  curl  u  =  0  (8.10) 

AT  =  0 

The  conventional  Hilbert  expansion 


e^'fn  (8.11) 


substituted  into  (8.2)  yields  the  sequence  of  compatibility 
relations : 


div  u  =  0 


VPq  ^   0 


(8.12) 


at  zero  order  and 


div  u  =0 
n 

Vp^  +  ijLq   curl   curl  u^_^   =  0  (8.13) 

AT^_1    =   0 


for  n  =  1,  2,  ....   We  recall  that  H-  =  n/e  is  a  scaled 
viscosity  coefficient  which  is  of  order  one. 


\?> 


Regrouping  the  variables. 


Pq 


Uq,  P^,  Tq  (8.14) 


gives  Pq  =  constant  to  lowest  order  and  a  well-posed 

Navier-Stokes  system  for  each  set  (u  ,  p   .,  ,  T  ) . 

n   ^n+1    n 

The  Stokes  linearization  gives  a  degenerate  result 
in  that  the  inhomogeneous  terms  vanish  at  any  order 
[this  is  easily  seen  by  induction  in  (8.8)].   Since  the 
same  differential  system  occurs  at  each  order,  the  sum 
of  any  finite  number  of  terms  in  the  Hilbert  expansion 
(8.11)  is  exactly  a  Navler-Stokes  flow;  e  enters  only 
through  the  boundary  conditions  which  must  be  found  out- 
side this  theory  and  which  we  expect  to  be  a  formal  power 
series  in  e. 

We  turn  to  the  Oseen  linearization  for  which  typical 
problems  are  the  flow  around  a  stationary  object  or  the 
flow  through  a  channel.   Taking  eU  as  the  unperturbed 


This  degeneracy  was  first  pointed  out  in  the  special 
case  of  one-dimensional  heat  flow  [Wang  Chang  and  Uhlenbeck, 
1955]. 
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flow  velocity,  the  formal  Hllbert  expansion  (8.11)  with 
the  same  grouping  of  variables  (8.l4)  gives 


p_  =  constant  (8.15) 


followed  by 


div  u_^  =  0 


U-Vu  +  Vp^  +  \i     curl  curl  u  =  0  (8.l6) 


I  U  .  VTq  -  AqATq  =  0 


then 


div  u-j^  =  U  •  VTq 
U-Vu^  +  Vp2+ |Xq  curl  curl  U-j^  =  a(U-V)VT  (8.17) 

I  U  •  VT-^  -  AqAT^  =  U  •  Vp^  , 


[  cj  =  ^  |X  +  §■  (  p,  ~   y-]  )/^  ]  ^"^^  further  inhomogeneous  forms 
of  the  same  differential  system.   The  system  of  equations 
(8.l6)  is  the  standard  incompressible  Oseen  linearization 
of  the  Navier-Stokes  equations  (with  scaled  viscosity  and 
heat  conductivity). 

In  the  general  case.  Re  ^   0,  kinetic  theory  enters 
more  substantially  than  in  the  special  case  Re  =  0,  where 
only  the  boundary  conditions  are  modified.   The  reason. 
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qualitatively,  is  that  the  real  characteristics  U  •  V 
of  the  Oseen  linearization  carry  the  boundary  values 
into  the  flow. 
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9.   Large  Reynolds  Number 

With  a  finite  value  for  the  unperturbed  speed  U,  the 
Reynolds  number  becomes  large  as  l/e,  and  there  can  be  a 
double  boundary  layer.   There  is  no  advantage  to  lineariza- 
tion in  this  generality,  so  we  may  Just  as  well  consider  the 
full  nonlinear  flow.   A  formal  Hilbert  expansion  yields 
Euler  equations  and  inhomogeneous  Euler  equations  as  in  the 
initial  value  problem.   This  system  is  mathematically  well 
posed  provided  that  all  streamlines  intersect  a  boundary 
(or  go  to  infinity).   If  we  adopt  the  direct  Hilbert  expan- 
sion, we  must  look  for  asymptotic  boundary  conditions  which 
are  appropriate  to  the  Euler  equations,  the  boundary  condi- 
tions to  be  extracted  from  an  independent  analysis  of  the 
combined  viscous  and  kinetic  boundary  layers.   A  boundary 

condition  on  u   (reflecting  the  displaced  boundary)  will  be 

1/2 
a  power  series  in  e  ^  .   The  first  order  term  is  taken  from 

a  classical  viscous  boundary  layer  analysis.   The  second 
order  boundary  term  0(e)  is  a  combination  of  a  kinetic 
contribution  and  a  second  order  viscous  effect,  and  it  is 
to  be  used  in  an  inhomogeneous  Euler  equation.   In  partic- 
ular, a  second  order  fluid  dynamic  boundary  layer  analysis 
is  not  correct  without  the  kinetic  modification. 

Another  procedure  is  to  give  up  the  attempt  to  examine 
the  whole  flow  and  concentrate  instead  on  the  viscous 
boundary  layer  alone.   The  stretched  variables  in  conventional 


52 


boundary  layer  analysis  Involve  e  (since  \i  ^  e) .      This 
reordering  of  terms,  when  combined  with  the  Hllbert 
expansion^ again  gives  Navier-Stokes  (in  a  boundary  layer 
approximation)  and  Inhomogeneous  Navler-Stokes  equations 
[Trilling,  196^].   Any  such  higher  order  boundary  layer 
analysis  must  be  matched  to  an  appropriate  fluid  formula- 
tion in  the  interior,  e.g.,  inhomogeneous  Euler  or  Navler- 
Stokes  . 

The  most  general  procedure,  applicable  uniformly  to 
the  entire  fluid  domain.  Is  to  follow  the  ad  hoc  rule 
described  previously.   The  complete  (nonlinear  and  nonsteady) 
Chapman-Enskog  expansion.  Involving  differential  equations 
of  increasing  order,  is  asymptotically  equivalent  to  a 
sequence  of  Navier-Stokes  and  Inhomogeneous  Navier-Stokes 
equations  with  inhomogeneous  terms  of  increasing  order. 
Analysis  of  the  kinetic  boundary  layer  should  give  boundary 
conditions  for  the  Navler-Stokes  equations  as  a  formal  power 
series.   The  number  of  boundary  conditions  does  not  change 
with  the  order  of  approximation.   Since  e  occurs  nontrivially 
in  the  Navier-Stokes  equations,  their  solutions  will  not  be 
simple  power  series  in  e.   Nevertheless  this  procedure  does 
produce  a  sequence  of  fluid  approximations  which  is  asymp- 
totic to  the  solution  of  the  Boltzmann  equation  to  any  order 
6   in  a  region  which  includes  the  viscous  boundary  layer  and 
excludes  only  the  kinetic  boundary  layer. 
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It  Is  Important  to  note  that  there  is  an  essentially 
arbitrary  choice  of  fluid  procedures  available  for  approxi- 
mating to  the  interior  fluid  flow,  excluding  the  kinetic 
boundary  layer.   The  number  and  type  of  boundary  conditions 
that  must  be  supplied  by  the  kinetic  boundary  layer  analysis 
depends  on  the  particular  fluid  choice  that  is  made  (Euler, 
Navier-Stokes,  Burnett,  etc)-  In  other  words,  the  kinetic 
boundary  layer  analysis  cannot  be  made  universally  but  only 
in  conjunction  with  a  particular  treatment  of  the  interior. 
The  common  worry  that  Burnett  and  higher  order  Chapman-Enskog 
equations  presumably  require  an  increasing  number  of  boundary 
conditions  poses  no  real  difficulty.   There  is  no  compelling 
reason  to  choose  the  Chapman-Enskog  procedure  which  is  only 
one  of  a  number  of  valid  formal  series.   The  choice  we 
recommend,  a  sequence  of  Navier-Stokes  systems  with  inhomo- 
geneous  terms  of  increasing  order  and  similar  slip  boundary 
conditions,  has  the  overwhelming  advantage  of  including  the 
viscous  boundary  layer  (as  distinguished  from  the  Hilbert 
expansion)  and  of  utilizing  a  familiar  system  of  differ- 
ential equations  (as  opposed  to  the  Chapman-Enskog  expansion). 
For  reasons  that  will  appear  in  the  next  section,  the 
inhomogeneous  terms  in  the  differential  equations  and  the 
finite  slip  boundary  conditions  enter  nontrivially  at  the 
same  order. 
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10.   The  Boundary  Layer--Sllp  Flow 

The  initial  layer  Is  studied  by  Introducing  a 
stretched  time  scale,  x  =  t/e.  In  terms  of  which  the 
Boltzmann  equation 


^+L(f)  =  -ee-||  (10.1) 


exhibits  weak  coupling  of  spatial  gradients  with  a  purely 
time -dependent  problem.   Similarly,  the  boundary  layer 
can  be  studied  by  Introducing  y  =  x-,/e  as  the  scaled  dis- 
tance from  the  wall  and  x'  =  (xp,x^)  as  unsealed  coordinates 
within  a  parallel  surface  x,  =  constant  to  obtain  an 
approximately  one -dimensional  Boltzmann  equation 


^1  11+  L(f)  =  0(e)  (10.2) 


This  estimate  [ignoring  the  term  0(e)]  Is  sufficient  to 
calculate  first  order  slip  coefficients.  We  can  easily 
refine  this  estimate 


4^  11+  L(f)  =  -ee'.  ^+  0(e2)  (10.3) 


In  order  to  simplify  the  calculations  (in  particular 
keeping  all  tensor  moments  of  f  unchanged),  we  do  not 
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transform  i,    keeping  the  original  rectangular  coorainates 
in  velocity  space. 

The  slip  boundary  conditions  arise  from  matching 
solutions  of  the  boundary  layer  equation  (10-3)  with  the 
Hllbert  expansion  In  the  interior.   To  do  this  we  shall 
have  to  postulate  certain  existence  theorems,  both  for 
the  full  Boltzmann  equation  in  a  general  Interior  domain  D, 
and  for  the  one -dimensional  equation  in  a  half  space.   The 
two  types  of  existence  theorems  are  quite  distinct.   For 
example,  u  =  0  can  be  imposed  on  the  entire  boundary  of  D, 
but  conservation  of  mass  implies  u^  =  0  at  one  end  (or  at 
infinity)  if  it  holds  at  the  opposite  end  of  a  one- 
dimensional  problem.   Also,  suitable  regularity  and 
boundedness  conditions  at  infinity  are  needed  for  the 
unbounded  domain. 

The  basic  boundary  condition  involves  specifying  the 
half  distribution  function  f~,  ^  •  n  <  0,  where  n  is  the 
outward  normal  to  D.   To  be  able  to  guarantee  mass  conser- 
vation (u  =  0)  at  the  boundary,  we  specify  f~  modulo  a 
constant  multiple  of  t/^q  [cf.  (6.6)] 

f-  ^  j6~  +  c^if^  (10. '^) 


The  parameter  c^^  (which  is  a  point  function  over  the 
boundary  of  the  domain)  is  to  be  adjusted  to  satisfy 
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u  =  0.   A  special  case  of  this  boundary  condition  Is 
diffuse  reflection  for  whj.ch  P~   Is  explicit, 

f-  =  CqI^q  +  u;  .  ^'  +  (3/2)^/2  T/i^        (10.5) 

Here  u  and  T  are  the  tangential  wall  velocity  and  the 
wall  temperature,  and  i/^  is  the  tangential  component  of 
the  vector  i/ ,    (6.6  ).   The  parameters  u'  and  T  are 

1  ^      '  ^  WW 

given  functions  on  the  boundary,  and  Cq  is  to  be  deter- 
mined such  that  Uj^  =  0.   This  is  the  linearized  version 
of  the  statement  that  the  emitted  distribution  is  a 
Maxwelllan  at  the  local  wall  velocity  and  temperature. 

We  now  list  the  conjectured  theorems  for  an  Interior 
domain  D: 

Theorem  1:   A  unique  steady  solution  exists  subject  to  f~ 
given  on  the  boundary  of  D. 

Theorem  2:   A  steady  solution  exists  with  u  =  0  subject 
to  (10.4)  where  j25~  is  given  and  Cq  is  left  open;  the  solu- 
tion is  unique  if  the  total  mass  is  specified,  /  p dx  = 

D 
given  (or  /  p dx  =  given). 

Corollary:   A  unique  solution  exists  subject  to  the  diffuse 
boundary  condition  (10. 5)  with  specified  wall  velocity  and 
temperature,  provided  that  the  total  mass  of  the  system 
Is  specified. 
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The  result  quoted  In  Theorem  1  Is  a  proved  theorem 
In  a  certain  weighted  Lp  norm.    In  the  same  norm. 
Theorem  2  and  Its  corollary  can  probably  be  proved  using 
the  methods  of  [Grad,  1966]. 

In  a  one -dimensional  steady  flow  there  are  five  con- 
stants of  the  motion.   In  a  nonlinear  version  these  are 


P"x 


2 


P"x  +  P  +  Pxx 


pu  u  +  p  (10.6) 

^^  X  y   ^xy  ^     ' 


pu  u  +  p 
•^  X  2   ^xz 


pu^(e  +  ^  u  )  +  u^(p  +  p^)  +  q^ 


and  In  the  linearized  version  (about  zero  free  stream) 
they  can  be  taken  as 


^x 


P  +  Pxx 


p    )      p 
^xy    *^xz 


^x 


(10.7) 


C  Cerclgnanl,  private  communication. 
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The  basic  conjectured  existence  theorem  Is 

Theorem  3:   A  unique  steady  solution  exists  in  a  half 

space  X  >  0  if  at  X  =  0,  f~  is  given  modulo  CqV'q  (10.4), 

and  if,  in  addition,  the  five  constants  (10.7)  are  specified. 

Far  from  the  imposed  boundary  condition  at  x  =  0,  we 
expect  a  Navier-Stokes  flow  with  linear  temperature  gradient, 
dT/dx  =  const.,  and  velocity  gradient,  du'/dx  =  const., 
u'  =  (u  ,u  ).   Setting  three  of  the  constants  equal  to  zero, 

y  2 


p    =  p    =  q   =  0  (10.8) 

^xy   ^xz   ^x  ^    ' 


guarantees  that  the  flow  is  bounded  at  infinity.   The  flow 
at  infinity  will  then  be  a  constant  Maxwellian  characterized 
by  five  parameters.   To  produce  a  flow  which  decays  to  zero 
at  infinity,  we  must  impose  five  additional  conditions. 
There  are  two  constants  (10. 7)  left  which  we  can  set  equal 
to  zero,  u  =  p  =  0.   This  leaves  three  finite  constants, 
u'  and  T,  at  infinity.   Subtracting  these  from  the  given 
boundary  function  f",  we  are  led  to 

Theorem  4:   There  exists  a  unique  steady  solution  in  the 
half  space,  vanishing  at  infinity  and  satisfying  the  boundary 
condition 

f~=p~+   CqI^q  -  u'--^'  -  (3/2)^/^  'Yii^  (10.9) 
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where  ^~   Is  arbitrarily  given  and  (cq,-u',T)  are  adjustable 
parameters. 

The  Interpretation  of  Theorem  4  Is  that  to  any  given 
emitted  distribution  j25"  are  assigned  unique  values  of  the 
slip  velocity  u'  and  temperature  T.   The  parameter  count 
(which  would  give  four  boundary  conditions  at  Infinity  for 
the  elementary  boundary  condition,  f~  =  given)  can  be  easily 
discovered  as  a  plausible  consequence  of  the  basic  Theorem  1 
applied  to  a  finite  one -dimensional  domain,  say  0  <  x  <  x-j^, 
on  letting  x,  — >  «.   The  application  of  Theorem  2  to  this 
domain  Is  somewhat  different.   Two  adjustable  constants  c^ 
are  available  at  the  two  ends  even  though  u  =  0  at  one  end 
Implies  u  =  0  at  the  other.  The  counting  Is  correct, 
however,  since  the  second  constant  Cq  Is  needed  to  adj'ust 
the  total  mass,  J  p  dx. 

We  shall  also  refer  to  the 
Corollary:  Theorem  4  holds  If  the  Boltzmann  equation  Is 
replaced  by  an  Inhomogeneous  equation 

4l  If  +  L(f)  =  il  (10.10) 

where  "^   decays  rapidly  as  y  — >  oo.   in  this  case  the  moments 
(10.7)  are  not  constant  but  are  taken  to  be  zero  at  y  =  00. 


In  a  special  case  In  which  the  parameter  count  Is 
simpler,  this  theorem  has  been  proved  for  a  model  of  the 
Boltzmann  equation;  Y.  P.  Pao,  private  communication. 
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The  slip  coefficients,  u'  and  T  In  (10.9),  will 
depend  linearly  on  the  given  boundary  function  ff   and  the 
Inhomogeneous  term  ii . 

Let  us  now  consider  a  general  domain  D  with  smooth 
boundary,  taken  with  the  diffuse  reflection  boundary  condi- 
tion (10.5),  u  and  T  being  specified  on  the  boundary; 
(the  boundary  may  be  Interior  or  exterior  and  may  have 
disjoint  components,  as,  for  example.  In  a  flow  between 
two  cylinders  or  spheres).   The  Boltzmann  equation  Is  to 
be  solved  In  the  original  scaling. 

As  In  the  Initial  layer  analysis,  we  Iterate  back  and  forth 

between  the  Interior  and  boundary  layer.   To  start,  we  solve 

the  Navier-Stokes  equations  in  the  given  domain  D  subject 

to  u  =0  and  the  given  boundary  conditions  for  u'  and  T  . 
n  D  .7  WW 

Using  this  solution  we  construct  a  Chapman-Enskog  function 
fg  +  ef-,;   fpi  is  a  local  Maxwellian  with  the  Navier-Stokes 
parameters,  and  f,  is  a  more  complicated  function  of  ^ 
(obtained  by  solving  Integral  equations)  involving  gradients 
of  the  Navier-Stokes  velocity  and  temperature.   The 
Maxwellian  function  f^  satisfies  the  kinetic  boundary  condi- 
tion (10.5)  exactly  (and  determines  a  specific  value  for 
the  parameter  c^).   The  difference  between  the  exact  solution 
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to  the  Boltzmann  equation  and  the  Navier-Stokes  approxi- 
mation, 

f  -  (fQ  +  ef^)  =  eg  ,  (10.12) 

satisfies  an  Inhomogeneous  Boltzmann  equation  (related  to 
the  Burnett  approximation)  and  Is  subject  to  a  boundary 
condition  Induced  by  the  Navier-Stokes  residual  f-,  , 

g"  =  c^l^Q  -  f^  (10.13) 

Turning  to  the  one -dimensional  boundary  layer  equation 
(which  to  this  order  is  simply  the  homogeneous  one- 
dlmenslonal  Boltzmann  equation),  we  find  a  solution  g„ 
(Theorem  4)  which  decays  at  infinity  and  satisfies  the 
boundary  condition 

g-  =  -f-  +  c^^Q  -  u[.if'    -    (3/2)^/2  T^1^4     (10.14) 

Having  determined  the  slip  coefficients  uj  and  T-,  (also  c-,), 
we  then  compute  a  Navier-Stokes  function  g-,  in  D  which  takes 
the  boundary  values  u-!,  T-,  (and  u  =0).   The  sum, 
gQ  +  g-|^  =  g  satisfies  the  boundary  condition  (10.13) 
(including  the  determination  of  c,  )  and  leaves  a  residue 
which  is  O(e^) . 
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The  first -order  slip  coefficients  u^  and  T-,  are 
computed  as  the  solution  of  a  one -dimensional  problem 
(Theorem  4)  with  a  boundary  condition  fJ  which  Is  taken 
from  the  first -order  Chapman-Enskog  function.   This  func- 
tion, f-.,    Is  a  linear  combination  of  velocity  and  tempera- 
ture gradients,  Bu./^x.  and  ^/Bx   (or  p..  and  q.),  taken 
from  the  zero-order  Navler-Stokes  solution.   A  scalar  T-, 
and  tangential  vector  u-I  can  be  combined  with  the  normal 
n.  as  linear  functions  of  p .  .  and  q^^  only  In  the  combina- 
tions 

T-,  =  ax  +  Pq 
1     n   ^^n 

(10.15) 
U-,  =  YT  +  6q 


where  t.  =  p .  .n .  Is  the  vector  force  on  the  wall,  and 

T  =T'n,  T'=T-nT  are  Its  normal  and  tangential 
n      '  n  ^ 

components.   The  coefficient  p  Is  the  conventional  tem- 
perature slip  due  to  normal  heat  flow,  and  y   describes 
the  slip  velocity  uJ  Induced  by  the  tangential  stress  t'. 
The  slip  u-I  caused  by  q'  Is  called  thermal  creep. 

The  teriH  ax  Is  ostensibly  of  higher  order  since  t 
always  vanishes  at  the  boundary  In  any  Navler-Stokes  steady 
flow  with  no-sllp  boundary  conditions.   But  It  would  be 


This  was  pointed  out  to  the  author  by  Y.  Sone. 
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retained  In  a  •clme -dependent  boundary  value  problem,  or 
In  any  problem  (such  as  with  the  thirteen  moment  equations) 
In  which  no  strict  ordering  with  respect  to  a  parameter  Is 
followed. 

The  formulation  for  second  and  higher  order  slip 
boundary  conditions  Is  quite  similar.   As  In  the  analysis 
of  the  Initial  layer,  we  note  that  both  the  Interior  and 
boundary  layer  expansions  Individually  satisfy  the  Boltzmann 
equation  without  Interference.   The  two  expansions  are 
coupled  only  through  boundary  conditions.   The  Interior 
(Hllbert)  expansion  can  be  written 


ftr  =  f Q  +  ef 2  +  e  f 2  + 


+  efQ  +  e^f]_   +  .  .  .  (10.16) 

+  e  f Q  +  ... 


where  f^,  f^,  f^  are  locally  Maxwell  Ian;  f-,  and  fp  are  the 
Navler-Stokes  and  Burnett  residues  Induced  by  f„;  f.  Is 
the  Navler-Stokes  residue  Induced  by  f„.   All  the  terms  on 
a  given  line  In  (10.16)  are  uniquely  determined  once  the 
leading  (Maxwell Ian)  term  Is  specified.   The  Maxwelllan  f^ 
Is  a  solution  of  the  Navler-Stokes  equations  subject  to 
the  given  boundary  conditions,  u'  and  T  .   The  higher  order 

WW 

Maxwelllans  f^  and  f^  satisfy  the  Navler-Stokes  equations. 
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but  with  slip  boundary  conditions  u-,,  T-,  and  Up,  T^  which 
are  to  be  found  from  the  boundary  layer  analysis. 

The  determination  of  f^  (from  the  first  order  boundary 
conditions,  u|  and  T-,)  Is  described  above.   The  second 
order  boundary  conditions  for  f^   are  found  (Theorem  4, 
Corollary)  by  solving  the  Inhomogeneous  one -dimensional 
boundary  layer  equation 


^^|fi+  L(gl)  =   -V  -^  (10.17) 

where  g  Is  the  already  known  first  order  boundary  layer 
solution  [called  g^   In  (10.14)],  and  g   Is  to  be  found 
subject  to  the  boundary  condition 

(g^)"--(f2+^i)"+^^0-^2-^'  -(^/^)^\-^l^         (10.18) 

The  Navler-Stokes  term  (f-i)~  In  the  boundary  condition 
leads  to  exactly  the  same  slip  coefficients  (a,p,7,B)  as 
In  (10.15)-  The  Burnett  function  fp  Involves  second 
derivatives  of  u  and  T   (or  first  derivatives  of  q.  and 
p.  .)  and  gives  rise  to  seven  slip  coefficients: 

T^  =  ak  +  ba  +  c(n.Vp) 
2     n     n    ^   ^'^ 

(10.19) 
u'  =  dk'  +  ea'  +  fV'p  +  gS' 
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where 

^Plr 


Sp„_ 

0  =  n  n  ^s-^-^  (a     =  e    )  (10.20) 

1    r  s  dx.    ^  n    n'      ^     ^ 

■^1  =  "r  ^ 


The  terms  In  Vp  are  alternative  forms  equivalent  to 

^p.  ./bx..     The  function  g   (generated  from  f-,  as  a  boundary 


condition)  takes  the  form 


(10.21) 


In  evaluating  ^'-Sg  /Sx'  for  (10.1?)  we  have  contributions 
from  4'*Sp.  ./Sx'  and  |'-Sq./Bx'  which  are  easily  seen  to  be 
included  within  the  format  (10.19)  (In  fact,  only  five  of 
the  possible  seven  slip  coefficients  arise  from  this  source, 
since  a  and  Vp  occur  only  through  the  combination  a  +  Vp). 
But  there  are  also  contributions  from  i^  •'djd/'dx^    and 
^'  •  ^t/'/Sx  which  are  derivatives  holding  i   fixed  in  a  given 

rectangular  coordinate  system,  whereas  the  components 

2     2  2 

i-,   =   ^'n  and  ^  +  ^  =  |^  X  n|   vary  with  the  normal  n  to 

the  boundary.   A  simple  computation  gives  the  result 
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«'-!!7=  (^1  -MsX^Aj^j)       (10-22) 


and  similarly  for  if,   where  )6-.    and  jd^   are  the  derivatives 
of  j6   with  respect  to  its  arguments  i-.    and  4p  +  4^ 
respectively  and 


K.  .  =  K'a.a.  +  K"b.b.  (10.2;5) 

describes  the  curvature  of  the  boundary  of  the  domain; 

k'  and  k"  are  the  principal  curvatures  along  the  principal 

directions  a.  and  b.  and  k.  .t .  is  the  geodesic  curvature  in 

-L  -I-  -^  d       u 

a  tangent  direction  t . ;  also  k'  and  k"  are  positive  if  the 
geodesic  curvature  is  toward  the  interior  of  the  domain. 
From  the  products  Pj_^K^g  and  0.^^^^   Iri  ^'  •  Sg  /Sx'  we  obtain 
seven  further  slip  coefficients, 

2    2  n    '^ d^n  2^rs  rs 

(10.24) 


I  - 


Ug  =  Yg'^  *^  +  62*^   "^  1^2*^"^  ■•"  ^2  ^'"^ 


where 


K  =  ^  K_  =  i(K'  +  K")  (10.25) 


rr 


is  the  mean  curvature  (note  that  k-t  =  k. .t.  and  K-q  are 
necessarily  tangential  to  the  boundary. )  The  four  coef- 
ficients Oo  J  Po*  "Vp'  ^2  ^^^    similar  in  effect  to  a,  p,  y,  5, 
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but  with  an  extra  factor  k. 

Since  the  second  order  boundary  layer  equation  (10.17) 

Is  inhomogeneous,  the  moments  (10.7)  which  are  zero  at 

y  =  00  are  finite  at  y  =  0  (cf.  Corollary  to  Theorem  4). 

The  normal  velocity,  u  ,  Is  therefore  different  from  zero 

[the  adjustable  constant  c  In  (10.18)  Is  required  to  make 

u  =  0  at  V  =  ool.   We  have  an  effective  displacement  of 
n        "^ 

the  boundary  given  by  the  Navler-Stokes  slip  condition 


"n  =  ^(^n  +  2KqJ  +  t>(a^  +  2KT^+||--K^3P^3)     (10.26) 


The  presence  of  only  two  slip  coefficients  Is  a  consequence 
of  dlv  u  =  0  which  allows  u  to  be  evaluated  In  terms  of 

^  ,00 

the  first  order  mean  slip  velocity  u'  =  I   u'(y)dy  [obtained 

ro 
g°(y)dy  In  (10.21)]  using  (b  u'-dx  for  paths  on  the 
u  ■' 

boundary.   In  computing  drag,  heat  transfer  at  the  boundary, 

etc.,  the  jump  In  the  values  of  other  moments  (10.17)  across 

the  boundary  layer  would  be  needed. 

The  slip  coefficients  g  In  (10.19),  ol^   and  "K^   In  (10.24)^ 

and  the  terms  2kt   -  k  p   In  (10.26)  can  be  considered 

n    rs  rs     ^ 

to  be  of  higher  order  just  as  was  pointed  out  for  the  first 
order  coefficient  a. 

In  addition  to  the  explicit  appearance  of  boundary 
curvature  in  (10.24),  it  Is  also  implicit  in  the  Navler- 
Stokes  expression  for  p. .  in  terms  of  Su./^x.  and  also  in 
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the  global  relation  between  T  and  the  boundary  values  of 

w 

q  and  k  through  the  solution  of  an  elliptic  equation. 

A  rough  estimate  can  be  made  of  the  first  order  slip 
coefficients  following  Maxwell;  for  diffuse  reflection 
(cf.  [Grad,  1958]) 

a  =   -  1/^  P  =  iir/&)-^^^^    0.627 

(10.27) 
Y  =  (7r/2)-^/^-^  1.235     6  =  -  1/5 

More  accurate  values  have  been  computed  for  p  and  7  using 
half  range  polynomial  expansions;  for  Maxwellian  molecules 
[Ziering,  I96O] 

p  ^   0.766 

(10.28) 
7  ^  1A8 

For  a  relaxation  model  of  the  Boltzmann  equation  essentially 
precise  values  have  been  calculated  for  p  [Welander,  195^; 
Sone,  1966a],  for  7  [Willis,  196O;  Albertoni,  et  al.,  I963], 
and  for  5  [Sone,  1966b], 

P  -^  0.757 

(10.29) 

7  ^  1.437  5  '^  -0.306 
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and  for  five  of  the  sixteen  second  order  coefficients; 
e-f  [Cerclgnanl  and  Sernaglotto,  1966],  (2b -c)  and  b 
[Tamada,  1967],  (72  +  2e  -  4f )  and  5^  [Sone  and  Yamamato, 
private  communication]. 


O.189A  (10.30) 


2b 

-  c  ~ 

-0, 

.1980?. 

e-f  ' 

'^2 

+  2e  - 

■4f 

^  2.55>^ 

62^ 

b  '^ 

0, 

.2927A 

The  normalization  (10.15)  In  terms  of  stress  and  heat  flow 
Instead  of  velocity  and  temperature  gradients  makes  the 
values  of  the  first  order  slip  coefficients  Independent  of 
the  essentially  arbitrary  parameter  which  enters  Into  a 
relaxation  model .   But  the  multiplicative  factor  for  the 
second  order  coefficients  depends  on  this  parameter,  i.e., 
on  the  definition  of  a  mean  free  path  [which  in  (10. 30)  is 

A  =  (8RT/7r)^/Vv]. 

We  have  found  by  consistent  formal  expansion  that  the 
Navier-Stokes  system  with  no-slip  boundary  conditions  is 
followed,  at  next  order,  by  an  inhomogeneous  Navier-Stokes 
system  (comparable  to  Burnett)  subject  to  the  slip  boundary 
conditions  (10. I5).   In  particular,  slip  boundary  conditions 
are  formally  Inconsistent  with  use  of  the  ordinary  Navier- 
Stokes  equations.   This  is  in  contrast  to  the  initial 
value  problem  where  the  first  occurrence  of  the  homogeneous 
Navier-Stokes  equations  is  accompanied  by  nontrlvlal  slip 
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Initial  conditions.   The  reason  for  the  discrepancy  is 
that  the  standard  diffuse  reflection  boundary  condition 
states  that  f"  Is  a  local  Maxwellian.   If  the  initial  func- 
tion f(C,x,0)  is  taken  to  be  locally  Maxwellian,  the  initial 
slip  also  recedes  one  notch,  and  it  appears  for  the  first 
time  at  the  Burnett  (or  inhomogeneous  Navier-Stokes)  level. 

We  have  presented  the  boundary  layer  theory  for  a 
Stokes  linearization  of  the  fluid  equations.   It  is  applic- 
able after  trivial  modifications  [viz.  restoring  ^p.  ./Sx. 
for  ^p/Sx.  in  (10.19)]  for  an  Oseen  linearization  and,  to 
first  order,  even  for  the  nonlinear  steady  Navier-Stokes 

equations.   In  the  nonlinear  case,  the  expansion  becomes 

1/2 
one  In  e  '       rather  than  e.   The  linear  first  order  slip 

theory  is  valid  since  the  kinetic  layer  is  thinner  than 

1/2 
the  viscous  layer  by  a  factor  0(e  ^  )  (we  require  that 

1/2 
Me  '   be  small).  The  entire  analysis  is  for  smooth  boun- 
daries; the  behavior  near  a  corner,  leading  edge,  etc., 
requires  separate  investigation.   In  a  time -dependent 
problem,  the  first  order  slip  boundary  condition  (10.15) 
is  unchanged  (after  leaving  the  initial  layer).   The  second 
order  theory  is  modified  by  the  occurrence  of  an  additional 
inhomogeneous  term,  -Sg  /^t,  in  (10.17)- 

Implicit  in  the  analysis  made  above  is  a  conjectured 
existence  theorem  for  the  steady  Navier-Stokes  equations 
in  an  interior  domain  D  given  the  vector  velocity  and 
temperature  on  the  boundary  and  the  total  mass  in  D. 
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11.   Infinitely  Strong  Shock  Profile 

We  recall  that  for  a  weak  shock  the  Navler-Stokes 
profile  Is  explicit  and  universal  [Taylor,  1910;  Grad 
1952], 


(2f(x)  =  tanh  ex/A*  (11.1) 


Here  ^   represents  any  of  the  flow  variables,  density, 
pressure,  temperature,  or  velocity,  suitably  normalized 
to  vary  between  -1  and  +1;  e  Is  a  shock  strength  param- 
eter, 0  <  e  <  1,  e.g. 


Pi  "  Po 


(the  subscripts  0  and  1  refer  to  the  front-upstream  and 
back-downstream  respectively);   A   is  a  specific  mean 
free  path.   Expressed  in  terms  of  the  viscosity  coeffi- 
cient 


A*=2|3Qt  (11.5, 


where  m  =  pu  is  the  mass  flow  constant  and  a„  is  a 
niimerical  constant  which  is  exactly  unity  for  Maxwellian 
molecules  (iJ.  «  T)  and  differs  from  unity  by  less  than 
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five  percent  for  all  potentials  between  Maxwelllan  and  hard 

1/2 
spheres  { |i  ^  T  '^  ) . 

The  weak  shock  profile  can  also  be  found  from  the 

Boltzmann  equation  by  a  formal  expansion  In  e.    Scaling 

ex  =  y  yields  a  Boltzmann  equation 


e||=|Q(f,f)  (11.^) 


in  which  the  parameter  e  appears  exactly  as  in  the  Hilbert 
theory.   But,  as  in  the  steady  flow  linearization  (Sec.  8), 
the  fluid  compatibility  relations  are  degenerate,  and  we 
do  not  obtain  the  standard  Hilbert  or  Chapman-Enskog 
sequence  of  fluid  equations.   The  steady  flow  Euler  equa- 
tions in  one  dimension  can  be  completely  integrated.   The 
compatibility  relations  are  therefore  essentially  algebraic 
rather  than  differential  equations.   The  lowest  order 
compatibility  requirement  (corresponding  to  a  local 
Maxwelllan  f  )  yields  a  constant  state.   To  next  order, 
the  variational  Euler  equations  are  a  linear  homogeneous 
algebraic  system;  the  condition  for  a  nontrivial  solution 
(vanishing  determinant)  is  that  the  zero  order  constant 
state  be  sonic.   We  recognize  this  as  a  bifurcation  equation; 


P.N.  Hu,  private  communication 
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in  searching  for  a  weak  shock,  we  are  looking  for  a  non- 
unique  solution  which  branches  off  from  the  trivial 
constant  state.   To  next  order  the  classical  G.I.  Taylor 
solution  (11.1)  is  found.   This  expansion  in  e  is  not 
strictly  comparable  to  Navier-Stokes ,  Burnett,  etc. 
since  the  compatibility  relations  are  not  a  sequence  of 
fluid  dynamic  equations.   But  a  similar  singular  pertur- 
bation applied  to  the  Navler-Stokes  equations  will  yield 
the  same  leading  term;  applied  to  the  thirteen  moment 
equations  it  yields  two  correct  terms,  etc. 

For  finite  shock  strength  the  Navler-Stokes  profile  can 
be  solved  numerically  by  Integration  of  a  simple  plane  direc- 
tion field  [Gllbarg  and  Paoluccl,  1953]  or  by  expansion  in 
e  [Grad,  1952]  or  even  explicitly  for  the  special  value 
aQ  =  20/21  [Becker,  1922;  Thomas,  19^^]-  A  paradoxical  early 
result  was  that  the  shock  thickness  for  strong  shocks  could 
become  arbitrarily  small  or  arbitrarily  large  compared  to 
the  mean  free  path,  depending  on  the  temperature  dependence 
of  the  viscosity  coefficient.   This  difficulty  was  resolved 
[Grad,  1952]  by  the  simple  expedient  of  introducing  the 
downstream  mean  free  path,  'A-,  ,  as  the  reference  length 


P.N.  Hu,  private  communication 
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Instead  of  the  upstream  mean  free  path  A^.   Not  only  is 
the  limiting  thickness  finite  in  terms  of  A  ,  but  there 
is  a  definite  limiting  profile  for  each  of  the  flow 
variables.   The  importance  of  this  choice  of  scaling  lies 
In  the  fact  that  ^-lA^  is  unbounded  in  a  strong  shock 
(except  in  the  special  case  of  hard  spheres). 

We  now  suggest  that  for  the  Boltzmann  equation  there 
is  a  finite  limiting  profile  at  infinite  strength  for  the 
distribution  function  f(^,x)  when  the  correct  scaling  is 

Introduced  in  velocity  space,  viz.  at  the  downstream  sound- 

1/2 
speed,  c-j^  =  (5/3  RT.,  )  '    .      With  this  scaling,  the  upstream 

Maxwellian  reduces  to  a  5-functlon  ( c^^  -.0)  centered  at 
u„/c-,  =  4//5  ;  cf .  Fig.  1  where  the  upstream  Maxwellian  is 
indicated  by  constant  f  contours  centered  at  u-,/c-,  =  1//5"  • 
The  crucial  point  is  that  the  entire  upstream  distribution 
Is  at  a  definite  thermal  velocity  when  measured  on  the  down- 
stream scale.   Even  though  the  mean  free  path  is  strongly 
velocity  dependent,  the  mean  free  path  of  an  upstream 
particle  when  considered  to  impinge  upon  the  downstream 
distribution  is,  within  a  factor  of  order  unity,  equal  to 
the  mean  free  path  for  self  relaxation  of  the  upstream 
distribution. 
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Based  on  this  qualitative  remark,  we  introduce  the 
Ansatz 


f  =  fg  +  f  (11.5) 


where 


f^   =   Pe(x)5(?  -U-)  (11.6) 


5    F5^^^'-^^   -0 


and  f  is  taken  to  be  nonsingular  on  the  downstream  scale 
(it  may  be  unbounded  near  ^  =  u^,  but  it  is,  nevertheless, 
an  integrable  function  and  not  a  distribution). 
For  any  cutoff  potential  we  write 


Q(f,f)  =  Q-^(f,f)  -  vf  (11.7) 

where 


Q 


^(f,f)  =  [  f 'f'  V  dax3?^  (11.8) 

v(?,x)  =  J  f  1  V  dax3|-^  =  V  +  Vg         (11-9) 

in  a  standard  notation.   In  the  expansion  of  the  collision 
term  Q(f-  +  f,  fg  +  f )  we  indentlfy  Q(fg>  fg)  and  vfg  as 
5-functions  (a  collision  between  two  molecules  in  fg 
remains  within  f^.)  and  all  other  terms  as  legitimate 
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functions  on  the  scale  f .   Consequently,  the  Boltzmann 
equation  splits  Into 


and 


?  H  +  Vgf  -  2Q^(fg,f )  =  Q(f,f)  (11.11) 


where,  in  v{i,x)    we  set  ^  =  u-. 


Vq(x)  =  v(Uq,x)  (11.12) 


The  mass  conservation  moment  of  (11.10)  is 


""oTT^   V5  =  0  (11-15) 


We  have 


X 

Pg(x)  =  Pq  exp(-  ^  r  v^dx)  (11.14) 


and  therefore  fg  (11.6)  is  given  explicitly  in  terms  of 
the  "moment"  v   of  f .  Substituting  this  expression  for  fe 
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Into  (11.11),  we  oDraln  a  nonlinear  Integro-dirf erentlal 
equation  in  f  alone.   This  Boltzmann-like  equation, 
provided  that  it  has  a  solution  satisfying  the  boundary 
conditions 


f(-oo)  =  0 

2 


Pi  (^-u,)' 


(11.15) 


(27rRT^)- 


gives  the  solution  to  the  infinite  strength  shock  profile. 
No  shock  strength  parameter  remains,  and  the  downstream 
parameters  p,  ,  u,  ,  RT-,  ,  are  easily  scaled  out  of  the 
problem.   The  only  parameters  that  remain  are  those 
implicit  in  the  molecular  force  law.   The  equation  (11.11) 
for  f  is,  in  appearance,  at  least  as  complicated  as  the 
original  Boltzmann  equation.   But  it  has  the  great  advantage 
of  being  nonsingular,  or,  at  least,  less  singular  than  the 
Boltzmann  equation.   In  particular,  this  will  allow  us  to 
use  naive  approximations  for  f  which  would  be  worthless 
for  the  evaluation  of  T. 

The  simplest  approximation  to  any  distribution  f  is 
the  locally  Maxwellian;  in  one  dimension  this  involves 
three  moments.   For  the  total  distribution,  f,  this  is 
too  crude  to  give  a  profile  at  any  strength;  it  leads  at 
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most  to  the  jump  conditions  and  a  discontinuous  shock  front 

Applied  to  f,  this  approximation  becomes  valuable. 

1   2 
Taking  the  three  moments  (1,  ^-.,    p  ^  )  of  the  full 

Boltzmann  equation,  without  approximation,  yields  conser- 
vation of  mass,  momentum,  and  energy, 


i 

hi 


^^  f  d|  =  m 


f  d^  =  P  (11.16) 


e^  f  de  -  Q 


where  m,  P,  Q  are  constants.   For  f^  these  moments 
(x-dependent )  are 


^1  ^5^^  =  PS^O 


i 
i 

J  I  ^1^^  fg^^  =  |P5" 


?1  fgd^  =  pgu^  (11.17) 


For  f,  using  the  locally  Maxwellian  approximation,  we  have 
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f  ^f  fdi   =  pu^  +  p  (11.18) 

Tl  ^^^^d^  =  I  pu^  +  I  pu 

where  p(x),  u(x),  p(x)  are  the  standard  moments  of  f. 
The  constants  m,  P,  Q  can  be  evaluated  at  x  =  -<»  where 

^  =   ^5' 


^  =    Po^o 

P  =  PqU^  (11.19) 

n   1   3 
Q  =  2  Po^o 


Adding  (11.17)  to  (11.18),  we  obtain 


PgUQ  +  pu  =  PqUq 


P5UQ  +  pu^  +  p  =  PqU^  (11.20) 


Eliminating  p  and  p  from  these  equations  leaves  an  equation 
in  u. 
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(p5  -  PQ^^O^^^^  "  ^^^0  "*"  ^0^  =  °  (11-21) 


If  p^  Is  not  identically  equal  to  p  ,  then 


4u^  -  5uUq  +  Uq  =  0  (11.22) 


This  algebraic  equation  has  one  solution  u(x)  =  u^  which  we 
dismiss  as  being  incompatible  with  the  boundary  condition 

at  plus  infinity  and  another  solution 


u(x)  =  ^  Uq  =  u^  (11.23) 


Returning  to  (11.20),  we  also  find 


^(Pn  -  P.) 


^0 


P  =  IT  (Po  -  P^K 


(11.24) 


In  particular. 


'^  =  ^=TB%=  "^1  (11-25) 


Both  u  and  T  are  constant  throughout  the  flow,  leaving  p(x) 
as  the  only  nontrivial  moment.  This  is  the  entire  informa- 
tion to  be  gleaned  from  the  moment  equations,  but  there 
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remains  Eq.  (11.15)  for  pg.   The  moment  v_  of  f  Is  a 
function  of  temperature  and  Is  proportional  to  p, 
Vq  =  a(T)p.   Since  T  =  T^  is  constant,  we  can  write 

Vq  =  a^p  (11.26) 

where  a  is  a  constant,  and,  between  (II.13)  and  (11.24) 
obtain 

dp.     4a 


dx     .Q 


This   has    the   solution 


-  ITT  (Po  -  Pe^Pg  (11-2^) 


Ps   -      ^^^    -y  (11-28) 

°        1    +    e   "^ 


and 


Pi 

P   =   ^-TT7  (11.29) 

1   +  e  -^ 


(p]_  ==  ^PQ^'    ^1^° 


Pi    +  Po^  "^  c        ^ 

P   =   P   +   P5   =  ^^-7-  =   Po[|  +  §  tanh  y/2]         (11.30) 

1+e-^ 


-   82 


where 

y  =  x/^  (11.31) 


and 


The  quantity  v  =  a^pT  is  the  collision  frequency  of  a 
molecule  at  speed  u   in  a  Maxwellian  f , ;  or,  it  is  the 

collision  frequency  of  a  molecule  at  speed 

1/2 
Uq  -  u^  =  3u^  =  (JRT-j^)  ^   in  a  rest  frame  Maxwellian 

at  density  p.,  and  temperature  T,  . 

The  net  velocity  u  is  easily  obtained  from  pu  =  pp,u„. 


u  =  u  ^  ^  ^  ^  (11.33) 

°  4  +  e"^ 


The  total  moments  p  and  RT  are  obtained  by  integrating 

over  f  =  f  +  f  5., 
0 

1  +  I  e-y 

p  =  p 


^  (1  +  e-y)(l  +  I  e'^) 


(11.34) 


1  +  ^  e-y 

T  =  T 


1  (1  .^e-i-)^ 


-  8J  - 


The  T  profile  is  flatter  at  the  downstream  end  than  any 

-2v 
other  variable;  the  departure  from  T-,  is  e  ■^    instead  of 

e~^.   This  is  particularly  interesting  because  a  very 
slight  overshoot  of  T  is  sometimes  found  at  this  end 
[Salwen,  et  al,  l^Gk] .      Our  approximation  is  too  crude 
to  decide  this  question.   But  the  higher  order  departure 
of  T-,  -  T  indicates  that  this  effect  may  be  quite  sensi- 
tive to  the  exact  approximation  used.   Separate  parallel 
and  transverse  temperatures  are  sometimes  defined;  they 
are 


1  +  e-y 
Ml  + 1-  e  y)^ 

(11.35) 
^1 

1  +  ^  e  y 


The  greatest  anisotropy  is  at  the  upstream  end  where 
T„/T^  ^   4. 

The  fact  that  this  very  crude  approximation  yields  any 
answer  at  all  Indicates  that  a  moderate  Improvement,  e.g., 
equivalent  to  thirteen  moments  (or  a  modification  such  as 
that  of  [Butler  and  Anderson,  I967])  should  yield  quite 
accurate  results.   Investigation  of  the  shock  equation 
(11.11)  for  the  behavior  at  x  — >  ±00  and  ^  — >  00  might 
suggest  better  approximations.   Also,  if  f  is  unbounded 
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at  ^  =  Up^  [this  appears  likely  from  the  term  Q,  (f,f^) 
in  (11.11)]  a  bimodal  approximation  to  f  would  seem  to 
be  called  for.   For  f  this  would  involve  a  sum  of  three 
terms,  two  peaked  at  u^  with  spreads  of  T^  and  T-,  and 
one  at  u-,  ;  this  is  essentially  trimodal. 

Our  three-moment  approximation  is  closely  related  to 
the  Mott -Smith  bimodal  assumption.   But  the  specific  rela- 
tion (11.1^)  for  p  replaces  the  more  or  less  arbitrary 
moment  equation  which  is  needed  to  supplement  the  Mott- 
Smith  prescription;  also,  the  fact  that  u  and  T  are 
constant  arises  naturally  rather  than  by  fiat. 

Elaborations  of  the  Mott-Smlth  approximation  ( the 
most  accurate  is  probably  [Salwen,  et  al,  196^])  are 
found  to  be  consistent  with  our  statement  that  the  correct 
length  scale  for  a  strong  shock  is  the  downstream  mean 
free  path.   It  is  likely  that  any  of  these  ad  hoc  approxi- 
mations with  enough  adjustable  parameters  is  fairly 
accurate.   But  an  equivalent  amount  of  effort  is  likely 
to  produce  much  more  accurate  results  if  applied  to  the 
infinite  strength  equation  (11.11).  This  has  yet  to  be  done 

The  separation  (11. 5)  into  a  "molecular  beam"  fg 
plus  remainder  has  been  introduced,  independently,  by 
I.M.  Scholnick  and  D-L.  Turcotte  (private  communication). 
The  remainder,  f,  is  approximated  by  using  the  Navler- 
Stokes  equations.   This  use  of  the  Navier-Stokes  equations 
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can  only  be  considered  as  an  ad  hoc  procedure,  possibly 
more  accurate  than  our  local  Maxwelllan,  but  less  easily 
assessed  than,  say,  a  thirteen  moment  approximation  to  f . 
In  the  present  context  the  Navier-Stokes  approximation 
for  second  and  third  moments  of  f  has  no  intrinsic  signif- 
icance, as  it  does  in  a  small  mean  free  path  expansion. 

For  the  determination  of  f,  it  is  sufficient  to  know 
that  fc-  is  concentrated  near  ^  =  u^.   But  as  a  separate 

problem,  we  can  investigate  the  structure  of  the  distribu- 

1/2 
tion  fp-  on  the  small  scale  ^  -  u_  ^  ^  ^'^n  ^    '  ^°^   hard 

potentials  (harder  than  a  Maxwellian  fifth  power),  self 

collisions  between  f,.  and  itself  occur  at  a  much  slower 

0 

rate  than  fp.  with  f.   On  the  distance  scale  that  f^  decays 
to  zero,  we  can  set  Q(fg,fg)  =  0.   Equation  (11.10)  can 
be  solved 


^5  =  ^0  "^P^-4 


Vq  dx)  ,  (11.36) 


and  we  find  that  fc-  remains  Maxwellian  while  it  decays. 

The  fact  that  the  approximate  profile  (11.50)  for  an 
infinite  strength  shock  is  a  hyperbolic  tangent  just  like 
the  exact  weak  shock  profile,  allows  us  to  make  a  very 
simple  Interpolation  for  the  shock  thickness  that  should 
be  accurate  to  better  than  a  factor  two  for  all  strengths. 
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The  Navier-Stokes  weak  shock  formula  (11.1)  extrapolates 
to  an  Infinite  strength  shock  thickness  (e  =  1)  of  just 
A  ;  (we,  of  course,  use  the  downstream  temperature  for 
|j,(T)  to  evaluate  A  ).   The  equivalent  thickness  from  the 
present  theory  is  2A  (11,30).   For  hard  spheres  a  simple 
calculation  shows  that  ^  «,  O.98  A  .   in  other  words,  the 
infinite  strength  shock  thickness  estimated  from  (11.11) 
by  a  local  Maxwellian  approximation  is  just  twice  the 
weak  Navier-Stokes  thickness,  X*,  extrapolated  to  infinite 
strength. 

The  hyperbolic  tangent  profile  is  exact  for  a  weak 
shock  but  can  only  be  expected  to  be  a  rough  approximation 
for  strong  shocks.   However,  the  very  strong  asymmetry 
that  is  found  in  strong  shock  profiles  using  a  relaxation 
model  for  the  Boltzmann  equation  [Llepmann,  et  al . ,  I962] 

seems  to  be  specific  to  the  model  itself.   The  creation 

J'  • 
f  f -,  ,  is  replaced  by  a 

Maxwellian  source  at  the  local  temperature,  which  approxi- 
mates T„  ^  0  for  the  downstream  part  of  the  profile.   This 


The  Navier-Stokes  thickness  for  hard  spheres  calculated 
exactly  for  an  infinite  strength  shock  is  some  ^0^  smaller 
than  the  weak  shock  extrapolation  [Grad,  1952,  Pig.  l4]. 
Empirically,  it  is  more  accurate  (as  well  as  simpler)  to 
extrapolate  the  weak  shock  formula  than  to  calculate  the 
correct  Navier-Stokes  profile. 
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is  quite  different  from  our  scaling  which  states  that 
molecules  are  created  with  a  velocity  spread  comparable 
to  T-,  for  the  whole  profile.   More  precisely,  the  correct 
upstream  limiting  values  of  the  moments  u  and  T  of  the 
residual  distribution  f  are  not  known  (they  turn  out  to 
be  u,  and  T^  in  our  locally  Maxwelllan  approximation  to 
f ).   Only  if  it  were  to  turn  out  that  u  — >  u„  and  T  — >  0 
for  the  actual  solution  at  minus  infinity  would  our 
Ansatz  (11.5)  be  qualitatively  comparable  with  a  model 
equation. 
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12.   Conclusion 

Taking  fluid  dynamics  as  our  model,  we  can  look  at 
the  theory  of  the  Boltzmann  equation  as  a  subject  which 
Is  several  orders  of  magnitude  more  extensive,  containing 
an  endless  variety  of  boundary  layers  and  non-unlfonn 
limits.   Enough  evidence  has  been  accumulated  to  indicate 
that  a  large  part  of  the  subject  Is  accessible  to  standard 
methods,  particularly  asymptotics  and  boundary  layer 
analysis.   But  some  aspects,  possibly  Including  noncutoff 
potentials,  nonself-ad Joint  spectra,  nonlinear  existence  in 
the  large,  may  require  significantly  new  approaches. 

In  trying  to  project  from  the  past,  we  can  distinguish 
three  qualitatively  different  eras:   transport  coefficients, 
ad  hoc  polynomial  and  moment  methods,  and  the  blossoming  of 
more  precise  mathematical  Investigations. 

Until  19^9  the  subject  of  kinetic  theory  was  essentially 
coextensive  with  the  Chapman-Enskog  theory.   The  ralson 
d 'etre  of  the  Boltzmann  equation  was  to  compute  transport 
coefficients  in  order  to  complete  the  classical  equations 
of  continuum  mechanics.   The  phrase  "solve  the  Boltzmann 
equation"  is  synonymous,  in  this  era,  with  the  operation  of 
solving  the  Hilbert  Integral  equation.   The  concept  of  the 
Boltzmann  equation  as  an  entity  in  Itself,  to  be  solved 
subject  to  initial  and  boundary  conditions  was  missing 
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(except  for  Carleman,  whose  singular  relation  to  the 
subject  was  only  revealed  posthumously  [Carleman,  1957])- 

The  break  with  this  preoccupation  with  transport 
coefficients  was  made  by  the  introduction  of  the  concept 
of  Interpolating  different  levels  of  description  of  a 
gas  [Grad,  19^9]  between  the  (at  that  time)  almost  inacces- 
sible Boltzmann  description  via  a  distribution  function  and 
the  macroscopic  description  via  fluid  moments  through  the 
Chapman-Enskog  hierarchy.   Polynomial  and  moment  approxi- 
mations gave  an  ad  hoc  scheme  for  this  interpolation. 

In  a  few  specific  problems,  viz.  linearized  and  one- 
dimensional,  high  order  polynomial  expansions  were 
carried  out  [Wang  Chang  and  Uhlenbeck,  1952  and  1953; 
Mott -Smith,  1954;  Pekerls,  et  al.,  1962].   More  accurate 
(but  more  special  and  more  complicated)  methods  of  param- 
eter fitting  were  the  introduction  of  half  range  polynomial 
expansions  [Gross  and  Zlerlng,  1958;  Lees,  1959]  and 
bimodal  approximations  of  varying  degrees  of  complexity 
[Mott-Smith,  1951;  Salwen,  Grosch  and  Zlerlng,  1964]. 
Quite  general  schemes  have  been  Introduced  which  perturb 
by  polynomial  expansion  a  suitable  weight  function  chosen 
to  fit  the  problem  [Kahn  and  Mintzer,  19^5;  Butler  and 
Anderson,  I967] .   All  these  techniques  are  essentially 
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ad  hoc,  and  the  ratio  of  accuracy  to  effort  (convergence 
is  unknown)  is  a  sensitive  function  of  the  cleverness  of 
the  performer. 

Mathematical  precision  has  developed  more  recently 
in  two  directions:   in  general  estimates  of  convergence 
(i.e.,  in  existence  theory)  and  in  precise  estimates  of 
the  behavior  of  solutions  in  special  problems  and  singular 
limits . 

The  text  of  this  paper  serves  as  a  summary  of  the 
increased  precision  which  has  developed  in  special  prob- 
lems and  singular  limits. 

The  early  existence  theory  was  confined  to  spatially 
homogeneous  problems  [Carleman,  1935;  Wild,  1951; 
Morgenstern,  195^]  or  gave  estimates  which  became  intol- 
erable after  one  or  two  collision  times  [Carleman,  1957; 
Grad,  1958]  or  one  or  two  mean  free  paths  [Willis,  I96I]. 
More  recently  a  strong  existence  theory  has  been  developed 
for  the  linearized  initial  value  problem  [Grad,  1963b]. 
Weaker  but  very  promising  results  have  been  obtained  in 
boundary  value  problems  [Cerclgnanl,  private  ccmmunlca- 
tlon],  and  a  start  has  been  made  toward  a  nonlinear 
theory  [Grad,  19^5;  Pao,  private  communication]. 
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The  Judicious  application  of  ad  hoc  parameter  fitting 
directed  by  the  growing  supply  of  precise  Information 
about  singularities  offers  the  most  promising  hope  for 
accurate  practical  calculations.   We  have  already  pointed 
out  the  application  to  the  fundamental  solution  and  the 
strong  Shockwave.   There  are  endless  variations,  such  as 
to  extend  the  precise  result  for  high  frequency  sound 
[Grad,  I966]  to  lower  frequencies  where  It  may  overlap  the 
normal  mode  analysis  from  the  other  direction  [Wang  Chang 
and  Uhlenbeck,  1952;  Pekerls,  et  al.,  19^2;  Slrovlch  and 
Thurber,  I963];  or  to  modify  a  half  range  expansion  to 
Incorporate  the  velocity  boundary  layer  at  |  =  0. 

Extrapolating  towards  the  future,  we  submit  the  fol- 
lowing as  a  suggestion  with  regard  to  the  most  promising 
directions  and  most  urgent  needs: 

1)  stronger  existence  theorems  for  boundary  value 
problems  to  try  to  bring  the  level  up  to  what  Is  known 
for  the  Initial  value  problem; 

2)  further  exploitation  of  the  abundance  of  known 
singular  and  nonuniform  limits; 

3)  development  of  numerical  techniques  to  be  able 
to  make  use  of  the  many  explicit  answers  (in  near  free 
flow)  which  are  given  as  multidimensional  Integrals, 
and  to  serve  as  numerical  experiments  for  comparison 
with  theory; 
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4)   study  of  infinite  range  potentials  (which  is  the 
real  goal  of  all  kinetic  investigations),  either  to  dis- 
cover properties  of  the  collision  operator  or  to  solve 
specific  elementary  flows. 

It  is  safe  to  say  that  kinetic  theory  for  monatomic 
rarefied  gases  has  reached  a  stage  where  the  obstructions 
are  mathematical  techniques  rather  than  physical  principles 
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